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ABSTRACT 


Begii'ining with the equations for conservation of mass^ 
conservation of momentunij and conservation of energy for tiie 
inviscid. , I two-dimensional or axisyminetric adiabatic flow of 

i 

an ideal gas, similarity solutions have been found which give 

2 

the flow field to order-6 about power-law bodies in the 
hypersonic limit where 6 is a body slenderness 

parameter. Some years ago, the hypersonic small disturbance 
equations were used to obtain ”zeroth-order ” similarity 
solutions for flow over power-law bodies. The second-order 
solutions, which reflect the effects of the second-order 
terms in the equations, are obtained by applying the method 
of small perturbations in terms of the body slenderness par- 
ameter 6 to these zeroth-order solutions. The method is 
applied by writing each flow variable as the sum of a 
zeroth-order and a perturbation function, each multiplied by 
the axial variable raised to a power. When these expanded 
variables are substituted into the flovj equations, a 
zeroth-order set and a perturbation set of four first-order 
ordinary differential equations is obtained, and the axial 
variation drops out. These equations are integrated 
numerically from the shock, v/here the boundary conditions are 

known from the Rankine-Hugoniot relations, toward the body. 

2 

The order-6 solutions which are obtained are independ- 
ent of the slenderness parameter 6 and thus are universal 



< < 1 . 


in that the^y apply for all values of 6 for which 6 
Ho V/ ever, except w he n t ti e b oci y p o we r - 1 a w e >: p o n e n t J. s e q u a ]. t o 
unity, the velocity functions, which form part of the solu- 
tion, have singularities at the body surface. These 
singularities are an effect of the entropy layer caused by 
the nose bluntness. Since the singularities are not removed 
by any of several methods tried, the solutions can only be 
applied away from the body surface. (It is suggested for 
future work that the singularities probably could be removed 
by applying the method of matched asymptotic expansions.) 

In comparisons with, the exact solutions for inviscld 
riow over wedges and circular cones, the order-6 similarity 
results give excellent agreement .for 6 less than about 0.4, 

corresponding to wedge or cone angles up to about 20^. Over' 

2 

an even larger range, the order-6 surface pressure predic- 

■ 

tions are superior to the Newtonian pressure lavf. The 
2 

order-6 results are a significant improvement over the 
zeroth-order results for body angles greater than about 12^. 
In comparisons with experimental shock wave shapes and sur- 
face pressure distributions for 3/4-power axisymmetric 

2 

bodies, the order-6 similarity solutions, give good results, 
considering that Mach number and boundary layer displacement 
effects are not included in the theory. For body fineness 

p 

ratios near two, the effects of the order-6 terms are 
significant only very near the body nose, whereas for a 



? 

fineness ratio near unity the order ~ terir-s has a large 
effect over alniost the entire hody , 

The order- 6^ similarity solutions are developed for 
infinite Mach number, but the derivation shows that they are 

/ 

compatible with shock-strength perturbation solutions, 

which Introduce Mach number effects. Also, while all results 

obtained are for no flow through the body surface (as a 

boundary condition), the derivation indicates that small 

amounts of blovjing or suction through the wall could be 

easily accomodated, • Finally, it is noted that a correlation 

suggested by Hornung for the shock wave shape and body 

pressure distribution can be applied exactly to all of the 

2 

flovj variables in the order-6 similarity solution form. 

This finding suggests for future work a possible refinement 
of the present derivation, using the local body or shock 
wave slope as the small parameter. 
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CHAPTER I. INTRODUCTION 

4 


A great deal of research has gone into investigating 
solutions to the small disturbance equations for hypersonic 
flow. One area that has received particular attention is 

i 

that of self-similar solutions for power-law profile bodies. 
While the effects of shock wave strength have been investi- 
gated in connection with these solutions, there apparently 
have been no reported efforts to Investigate the effects of 
neglecting the second-order terms of the complete Inviscid 
flow equations in order to reduce them to the small distur- 
bance form. The purpose of the present study is to determine 
the effects of retaining these terms by using , a perturbation 
analysis to obtain second-order similarity solutions for 
power-law bodies. 

Since this dissertation will be concerned with finding 
a particular set of similarity solutions of the inviscid 
flow equations, it is important at the outset to establish 
what is meant by similarity solutions in hypersonic flow. 

The similar solutions referred to here are solutions for 
self-similar flows; i.e. flows for which the flow field 
(expressed in suitable coordinates) at any one position 
along the body is the same as that at every other position. 
(In the corresponding unsteady self-similar flows, the flow 
field in suitable coordinates at any one time is the same as 
that at every other time.) Inviscid axisymmetric supersonic 


1 



2 


flov/ over a cone with an attached shock v/ave is a classical 
example of a self- similar flow and represents a particular 
case of the similar solutions discussed herein. For the 
cone, the flow field properties (eg. the pressure, the density, 
and the velocity components) are themselves constant along 
rays from the cone vertex. For the other power-law bodies, 
the flow field properties are not constant themselves, but 
similarity functions describing these properties are 
constant (to the order of the solution) along curved power- 
law paths from the nose of the body. 

The similar solution approach to solving the flow 
equations is valuable because it allows a reduction in the 
number of independent variables in the problem. In particular, 
for hypersonic flow about power-law bodies, the similarity 
approach reduces a system of partial differential equations 
to a system of ordinary differential equations. As noted by 
Hayes and Probstein (1), generally these flows occur only 
for a self-similar fluid (the most practical example of which 
is a perfect gas with a constant ratio of specific heats) 
and a self-similar shock wave, i.e. one having the same 
density ratio across it at every position. 

All of the early investigations of similar solutions 
related to the present problem were concerned with unsteady 
flows. Early in the Second World War, Taylor (2) developed 
a similar solution for the flow behind the spherical shock 



wave produced by the instantaneous release of energy at a 
point (e.g.j an atomic explosion). Sakurai (3j^) generalized 
Taylor's approach to obtain solutions for cylindrical and 
planar shocks as well. He also introduced perturbation 
analysis as a means of obtaining solutions for more moderate 
shock wave strengths. The equivalence of these unsteady flows 
to steady flows in one additional space dimension was pointed 
out by Hayes (5). This equivalence applies to the inviscid 
flow equations reduced to the hypersonic small disturbance 
form, as derived by Van Dyke (6). 

Lees (7) found that there are self-similar flow fields 
for bodies having power-law profiles, and Lees and Kubota 
(8) determined the range of power-law exponents for which 
the similarity holds. Kubota (9) obtained numerical solutions 
for this case (herein called the "zeroth-order" case); he 
also applied a perturbation in the strong shock parameter (as 
Sakurai had done for unsteady flow) and numerically obtained 
first-order similar solutions for moderately strong shock 
waves. Mirels (10) computed additional and more accurate 
numerical results for the zeroth-order and moderately-strong 
shock wave cases. He also derived approximate analytical 
solutions for these cases. 

The parallel but independent work of Investigators in 
the USSR has been thoroughly described by Hayes and 
Probstein (1). Beginning at about the same time as Taylor, 
Sedov (11) studied the intense spherical explosion problem 



in a more general form and developed an analytic solution 
for it (12). Grodzovskii (13) and Chernyi (1^) applied the 
unsteady results to the steady hypersonic flow problem, 
Stanyukovich (15) and others Investigated a number of 
related problems. 

All of the important developments in the use of hyper- 
sonic small disturbance theory to obtain solutions for power- 
law bodies were treated in a unified way by Mirels (16), who 
added an analysis of perturbed power-law body shapes. More 
recently. Freeman (17) investigated the effects of the entropy 
layer caused by the nose bluntness of the power-law bodies 
and determined the power-law exponents below which the 
entropy-layer effects predominate. Again independently, 

Sychev (18) developed a correction to the power-law body 
shape to account for the effect of the entropy layer. 

A few experimental investigations of the flow field over 
power-law bodies have been made. Kubota (9) compared his 
theoretical results to surface pressure distribution and 
shock wave shapes measurements for 2/3” and 3/^”power bodies, 
obtaining good agreement for the more slender bodies. 

Peckham (19) measured pressure distributions and shock wave 
shapes for a series of power-law bodies, some of which fall 
in the similar-solution range. Freeman, Cash and Bedder (20) 
and Beavers (21) also presented detailed shock shape data 
for series of power-law bodies, registering some disagreement 



with Kubota's results. Spericer and Fox (22) present aero- 
dynamic drag and other data for several power-lav/ bodies 
over a wide Mach number range, Ashby (27) presents aero- 
dynamic data for a similar series of bodies over a range of 
Reynolds numbers at Mach 6, and Ashby and Harris (28) use 
method of characteristics and boundary layer computer pro- 
grams to show the Important effect of boundary layer transi- 
tion on the total drag of those bodies. 

Townsend (23) applied the zeroth-order solution of 
Kubota and Mirels^ with their shock-strength parameter per- 
turbation and a boundary layer displacement correction, to 
the problem of estimating the forces and moments on a half- 
axlsymmetric body under a thin, flat wing. In order to study 
a range of configurations at a moderately hypersonic Mach 
number, Townsend applied his method to configurations which 
are marginally slender, (i.e. to configurations for which j 
the errors arising from body thickness are small but not 
negligible). This type of application points up two 
reasons for seeking solutions which include the effects 
of the second-order terms for body slenderness in the flow 
equations: (1) to assess the error caused by making the 

small disturbance assumption, and (2) to improve the 
accuracy of calculations for marginally slender bodies. 
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When compared with experimental data for axisymmetric 
power-lavj bodies and for viirig — conical-body configurations, . 
Townsend's method gave good agreement where the basic 
assumptions were satisfied. An example series of computations 
with variations in the principal parameters at a full-scale 
flight condition showed that varying the power-law exponent 
has a greater effect on longitudinal stability and trim than 
on the lift-drag ratio. The computations for Mach 6 gave 
higher maximum lift-drag ratios, higher drag coefficients at 
zero lift, but essentially the same stability characteristics 
as their counterparts for Mach 12. 

In the present study the second-order similarity solu- 
tions were obtained by a perturbation method. This method 
used expansions of the variables in terms of a small parameter 
to obtain higher-order solutions as perturbations from a 
known zeroth-order solution. The approach was very similar 
to that of Sakural (4), Kubota (9), and Mirels (16) in their 
first-order determinations of the effects of shock wave 
strength; but, the small parameter used herein was a body 
slenderness parameter rather than the shock strength parameter. 
Van Dyke (.29) describes the application of perturbation 
methods to fluid mechanics, and Van Dyke (30) shows how, in 
favorable cases, such solutions can be extended to improve 
convergence when the perturbation quantity is not small. 



The importance of the results to be obtained from the 
present study lies in their practical application. The 
principle area for this is in estimating the aerodynamic 
characteristics of generalized configurations (e.g., 
Townsend’s (23) family of wing-body combinations). By 
improving the results of such studies and by better defining 
their limits of applicability, the present work contributes 
to their usefulness in suggesting designs (or parts of 
designs) for such hypersonic vehicles as transports or re- 
entry spacecraft. 

The remainder of this dissertation will describe the 
development of the solutions and present the results. 

Chapter II gives the theoretical development. It goes 
through the transformations of the flow equations required 
to put them into similarity form, discusses the results of 
keeping terms of different order, describes the application 
of the boundary conditions, and develops an alternative 
formulation of the problem. Chapter III presents the general 
scheme for solving the equations and deals with the diffi- 
culties which arise. Chapter IV discusses the results and 
their region of validity. Chapter V gives the conclusions 
reached as a result of this study. The Appendix describes 
an approximate analytical solution used near the body sur- 
face, where the equations are singular. 



CHAPTER II . THEORY 


A . Transformation of Basic Flow Equations 
This section will show how the basic flow equations can 
be transformed to obtain a separation of variables for the ' 
case of ' hypersonic flow over power-law bodies. The starting 
point for this process is the system of steady, two-dimen- 
sional, inviscid flow equations for a perfect gas in physical 
variables: 


Continuity : 

9pu 

3x 

3PV + PV = 0 

3r r 

Longitudinal 

momentum: 

u ^ 
3x 

+ v|^+i|f=0 
3r p 3x 

Lateral ■ 
momentum: 


+ vp+ip=0 

3r p 3r 

Energy : 

(“ h 

+ V “ 0 


The constant a in the continuity equation has the value 
0 for planar flow (Cartesian coordinates) or the value 1 
for axisymmetric flow (cylindrical coordinates). The bars 
over the variables indicate that they are dimensional 
quant Itites. 

Normalization . The initial treatment of these equations 
follows that of Kubota (9) (also covered by Mirels (l6)), 


8 
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except that no terms are dropped. Kubota showed that for 
slender bodies in hypersonic flow the variables can be norma- 
lized using the expressions: 


X 


X 

I ’ 





u = 


u - 

<S u„ 



) ( 2 ) 


The 6 is a body slenderness parameter (to be discussed 
later) Introduced so as to make the dimensionless variables 
of order unity. These variables are substituted into equa- 
tions (1) to obtain the normalized flow equations: 

^ 3pu . 9p . 3pv . ^ pv _ ^ 

Continuity: 6 


Longitudinal 

momentum; 


Lateral 
momentum : 


3v . 3v , 3v , 1 ^ 


) ( 3 ) 


Energy : 





Note that each of these equations contains a leading term in 

p 2 

6 . If the body were sufficiently slender, the order-6 terms 
could be dropped, leaving the hypersonic small disturbance 
equations used by other workers. For the present study, 
however, the equations are retained in the complete form. 
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Similarity variables . The next step is to put the flow 
variables into similarity forms. Still following Kubota (9), 
these will be found by comparison with the flow through an 
oblique shock wave. The normalized flow variables just 
behind an oblique shock are (2^): 


Po = 


.2- -2 
>5 


2 2 
Y6 M 


6 7M 


2 P, 


sin^Gg ~ (y-D 


Y + 1 




(y+1)M^ sln^G 
’ s 

Poo (y“ 1)M^ sin^G + 2 


u = 


u - u 
s « 


6 u 




6u. ^ '' 


,2(wf sin^Og - 1) 


-2(M^ sin^0„ - 1) 


(Y+1)M‘ 


(y+i)m‘ 


cot 0. 




If the shock wave shape is given by R(x), its slope is 

R»(x) = ^ ~ 6R'. The shock wave angle 0^ is 

related to slope by tan 0 = R’ = 6R^j from which 
2 ? ^ 


sin 0^ = 


Putting these results into the oblique 


2 2 ’ 

1 + 6 R' 

shock relations gives (for R* of order unity): 


^s Y+1 


R 


.2 


- Xrl 

l+6^R’^ 







la lia 
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Y + 1 
P = - — r 
s Y“1 


1 + ^ 


Y-l 


R 


,2 


= ^ + 0(6S 

(y-1) R' 


u_ =1- 


Y+1 


R 


.2 


2 2 
1+6"^R» 


-* e 


2 ^,2 ^ 2 . 2 „ t 4 


Y+1 


R 


Y-H 


v_ = 


s Y+1 


R 


+ 62r.2 R'J 


2 „ _2_ j2fj,3 2_ e_ . 


Y+1 


R’ - 


Y+1 


Y+1 R 


- + 0(6") 


(4b) 


‘(4c) 


(4d) 


Here e h ^ is a shock strength parameter; as e -> 0, 

£ -> the limiting value for shock wave strength. Using 

these equations as guides, the flow variables are taken to 
have the forms : 


p = ;ji^(r) + R‘^ + (r ) (R^”^ 

u = + 6^V2(r)R’^ + e v^(r) 




) (5) 


j 



At this point, in order to get an expression for the 
shock wave shape, consideration is narrowed to flows about 
pov/er-law bodies. Under the hypersonic small disturbance 
assumptions, a power-law body - x^) produces a power-law 


shock wave (R ~ x^) for 


3+d 


< m < 1. (See Lees and Kubota 


(8).) Specifically, for 6^ << 1 and e << 1, the ’’zeroth- 




order” shock shape about a body 


(1)"^ 


IT 


= 6 


= — (I) 

2f \ll 


is given by 


in normalized coordinates, the shock shape 


i • 1 m 

about a body r^ = 2 ^ ^ 


is = X 


ra 


(Figure 1, parts 


!(a) and (b)). Note that for m=l the body is- a wedge (for a = 0 

j ■ ■ . - . . . ^ . 

or a cone (for a = 1), both of which are known to have straight 
shock waves and therefore satisfy the above relations. For 

J ', ■■ ■' .. * ■ 

fn’'< “X, the power-law body has a small blunt nose, so that' the 

j 

shock wave is detached. Consequently, this type of relation 
between the body and the shock cannot hold in the immediate 

I 

vicinity of the nose. The effects of nose bluntness on the 

Iflow downstream are confined to a thin layer near the body 

1 

surface. Freeman (17) found that the effects are less than 
order for m > . (For y = 7/5 this amounts 


2^ 

to ^ ^ .77 for a = 0 and to 


2^ 

3^ 


. 63 for 0 = 1). 


: The expression for the zeroth-order shock wave shape 

serves to define the slenderness parameter 6. That is, 
-\m 


- - lx\^ 

Hq(x) = evaluated at x = £ gives the relation 

Thus, 6 is the tangent of an angle defined by 


j R i^) 

I o 

a = ° 


i. 


Thus 



Loterol coordinate, 77 Loterol coordinate 


Longitudinof distance, x 

(a) Physical coordinate system. 



Longitudinal coordinate, x 

(b) Normalized coordinate system. 



(c) Similarity coordinate system, 

Figure 1, Power-law and Zero-order Shock in Physical and 
! Transformed Coordijiate Systems, 
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the shock v/ave position (Figure l(a.)) and is, in fact, a 
'’mean shock v/ave angle” parameter. However, since the 
zeroth-order case implies (as will be seen) , the 

shock lies near the body so that the shock wave angle and 
body slenlerness are closely related. The definition of 6 

I 

is made in terms of the shock for convenience, since the 


solutions are to be found by integrating from the shock to 
the body. 

i 

I To aid in the separation of variables, a shock-oriented 

i 

coordinate system is now introduced (Figure 1(c)). This 
system has 


I 



i 


so that r = riR^ = The body surface is then 

^b “ ^ where = 2^* 

, The shock wave shape to be used in equations (5) is the 
zeroth-order shock with Kubota’s (9) shock strength pertur- 
bation and a separate perturbation for the body slenderness. 
It is taken to be 


! R(x) = R^(l + + ea^m^R^"^), 

i ■ 
i 

where the constants and a^ are to be determined as 

i 2 

part of the solution. (The factor m is included in the 



last term of this equation so that it conforms to the usage 


of Kubota and Mirels,) Putting in 


R(0 = C 


m 


1 + + sa,e2(l-«0 


( 7 ) 


j Substituting the derivative of equation (7) into equa^ 
tlons (5) and ordering the terms by powers of 6 and e 

p 

give expressions for the flow variables in the following form 
(neglecting terms of order 6^ and of order e^): 


(n)m^ ^ 2(l-m) ^ ^2 p ^ 




+ eF^(n >m' 


p(Cjn) = ^ 2(l“-m) 


u(^,n) = V^(n)m^ + 6^V2(n)m^ ^-^(1-m) 


+ ev^(n)m' 


(C,n) = <j>Q(n)m C K 




>( 8 ) 


J 


The relations between the functions of n and the functions 
of r in equations (5) are not needed since it is easier to 
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work with equations (8) directly. These expres.sions are now 
ready to be substituted into equations (3) to obtain the 
transformed flow equations. In taking the derivatives of 

these variables it is necessary to note that, whereas ^ - x 

. 3 f ^ -m 

and = Oj n = rx 

Bri _ -m _ r*~rn. 


{ 

gives 

i 

the 

simple 

relations 

3C 

3x 

= ; 

( 

i . 

gives 

1 

i 

i 

In 

9x 

= -mrx” 

m-1 

mn 

and 

i 


3 

= ii 

^ + 

9n' 

3 

j 


3x 

3x 


3x 

3n 

and 







1 

{ 


3 


^ + 

3n 

a 

i 

1 


3r 

3r 


3r 

3n 


3r 


= x “ = C ; thus 


) ( 9 ) 


L L- 

,m 3q 


I Equations in similarity variables . When the expressions 
i(8)”are substituted into the normalized flow equations (3)^ 

j 

i . . . 

these become (away from the nose ^ = 0): 


Continuity : 


( d) * + — d))~(ri-<j) 

^0 ^0 n o ' ' ^0 ^0 




:+ + ('f'o + ^ + [f’i + ^ ‘<’0 " 2 

\ ■ ■ ■ 

1 

- (n - - [nv- + 2 (^) vj 

;+ ejvi + (f; + 1 Vh * ^ + 2 4-1 


'(10a) 


- (1 - }i 


+ 0(6^) = Oi • 
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Longitudinal momentum : 


F’ ■ 

(n - + n + 2(^^)(v^ + -™) 

o o lii m o ;!; 

o o 


m3g-3+2m 


’ ( F ’ P 

+ 6^|(n “! + n + 7 ^) - 

I I ‘ ^ ^ m 2 0^2 


f ^[nF' + 2(~ 

I 0,2 o ' m 

t- 


)F. 


m o 


'^2 + 


nv' + 2(~)v 
o m o 


5--5+^m 


+ e5(n - *o)v' + n ^ ^ 

° 'I'o ^ 


nF‘ + 2(i-^)F 
o 'mo 


ilj. 


7- V 


^ + 0(6^) =' 0; 


Lateral momentum: 


f ^ 


(n - - 7^ + (~H 

FI F» 


2j_~2+m 

m 5 


6‘j(n - ^ ^ 

° '*'o 


(j)> _ 3(i-51) 
m 




+ (izSL)^ 




+ e 


j(n - ^ ^ (^3 


-m 


+ 0(6^) = 0; 


:(iob) 


(lOc) 


f 
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Energy : 




^_o 

> 

o 


F» 

o 


2 (^) 


o 


rn 


rn^' 


+6 (n - 


i|i' P' ilipt' 

o \|j 


p F . 

2 Q 


L pflzHL’i fv _ fv ^ 

i m ^ A - P„)4’2 




>p: 


p» 

o 


o J 


vjm3r3-^2m 


f e (n - 


(y 


'•'i 

♦I 


F{ ’('I't’o 

^) - (y 


F-.P’ 
1 o 


r 2( 


4^- 


P. 




^o o ^o 


F» 

P 


KlOd) 


+ 0(6^) = 0. 


•These equations are seen to be ordinary first-order differen- 

T ' • 

tial equations, linear in the derivatives of the functions 

! 

defining the pressure, density, and velocity fields. It is 
noteworthy that although 6 appears to the first power in 
the normalization of variables (equations (2)), only the even 

powers of 6 appear in the final form of the flow equations. 

i 

Thus, while the solutions to be found are of second order 

1 

i 

in the body slenderness parameter 6, they could be considered 

! ' 2 

of first order in 6 . To avoid any ambiguity they will 

j p 

generally be referred to as order~6 solutions. Physically, 

i 

the absence of lower powers of 6 indicates that the error 

due to a given body thickness is less at hypersonic speeds 

’ 'R/2 

than at lower speeds, where terms of order 6 or 6 

(Van Dyke (6) ) , 


appear 



/ 


B. Discussion of Orders of 


As a result of the normalization procedure (equation (2)) 
the variables p, p, u, and v are of order unity for slender 
bodies in hypersonic flov/. The similarity variables in 

I 

pquations (10) are also considered to be of order unity, but 

i ^ 

phis assumption must be tested by the results obtained. The 
pevelopment so far has been based on 6^ and e being small 

I 

parameters. This section will consider their relative 
sizes. 


i Zeroth-order equations . If 6^ << 1 and e << 1, so 
phat all terms containing either one may be neglected, 
equations (10) are reduced to the zeroth-order equations: 




bontinuity: ^ ^ ^ ^ ^ 

I F» . p 

Longitudinal (n - (|) )v’ + n -^ + 2(i^)Cv + 

momentum : 9 ^ ^ ^ 


= 0 


Lateral momentum: (n - rf> >rf)’ - — + = n 

r — '^•o^^o ^ m ^ 

f 

Energy : 


) ( 11 ) 


(n - ^ - 2 (i^) = 0 


^O o 


J 


These equations represent the case first studied by Kubota 
j(9)i they are the same as his equations except that he 
pmltted the longitudinal momentum equation, which is 
uncoupled from the others. References 9, 10, and 23 contain 
results of numerically solving Kubota’s equations, which are 
a special case of the present more general treatment. 



20 


Equations (11) contain only one parameter, the povjer-law 
exponent m. Thus, for two-dimensional flov/ (a = 0) or 
axisymmetric flow (a = 1) of a given gas, the similar solu- 
tions P^(rj), v^(n), and each form two 

families jof "universal functions" depending only on the 
i 

power-law of the body. 

As was mentioned in Chapter I, the simultaneous applica- 
tion of the conditions leading to these equations imposes a 

2 

stringent condition on the Mach number; viz., 6 << 1 and 

e = — << 1 requires M >> 1. This relation between c 


and 6 is illustrated in Figure 2, where part (a) shows, 

2 

for example, that e < ,1 and 6 < .1 are both true at 
Mach 12 only if 6 - .3* Note, in addition, that dropping 
the terms in e removes all Mach number dependence from the 
equations, which really implies » and Illustrates 

Hayes' "Mach number Independence principle" (ref. 1). 


Order-e equations . If, in equations (10), the terms in 

2 

6 << 1 are dropped but the terms in c are retained, two 

systems of equations can be obtained by setting the zeroth- 
order and order-e terras separately equal to zero. The 
zeroth-order system is the same as before; the order-e 
system is: 


Continuity : 


Ui 6 ' + 

^ 0^1 




n ^0 ^1 


(n - + [(j)^ + 



+ 2 ( 


1-m 

m 


(12a) 


0 








SHOCK STRENGTH PARAMETER 



(c) Relative errors from neglecting terms of order and of order 6^. 

Figure 2. Concluded,. 


r\3 












'■’I 




nB" + 2(i^)F 
o m o 


^1 - = 0 


Lateral momentum; 


|(n - ^ ^ ^ + (i^)] ^ = 0 

) O . - 


Energy ; (n - <j>^) 

1 V/ 


(y 




pt pi 

o iB P . 

o 


T ^ Ft F' 

r ^ - jl) - (V / - p^) h = ° 

o ^o o 


(12b) 


(12c) 


(I2d) 


These equations are the same as Kubota’s (9) first-order 

i ' ■ . . • • 

15'erturbati on for shock wave strength, except that (again) he 

I 

omitted .the longitudinal momentum equation since it is un- 

coupled from the rest. They can be solved numerically using 

1 

the results of the zeroth-order solutions. The resulting 
similarity functions F^(ti), Vj^(n)j and (|>^(n) are, 

i 

like the zeroth-order functions, universal in that they 
depend only on m as a parameter. References 9, 10, and 23 
contain the results of the numerical solution. Applying this 

shock wave strength perturbation reintroduces the Mach number 

i 

I 

dependence and somewhat relieves the requirement that • -> » 

j 

|(see Figure 2(b)), but the body must still be very slender ■ 

i . ’ 

■ 2 
in order that 6 << 1. 
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2 ' 2 
Order~6 equations . If the terms in 6 are kept in 

equations (10) and all higher order terms are dropped, a case 

is found which has not been studied previously. This is the 

case of present interest. Since general similarity solutions 

1 ' . ■ . 

are being sought which do not depend on the particular values 

i 

of 6 or e, each of the three major terms in each of the 

four conservation equations must be separately equal to zero. 

As a result of observing this, the terms can be separated into 
! ■ ' ' ' 
twelve equations in the twelve unknown functions F^, 

I 

■ ^1’ ^1" ^1’ ^1" ^2’ ^2’ ^2" ^2‘ of* these 

equations are the zeroth-order and order-e systems of equa- 

‘ 2 
tions found before. The remaining four are the order-6 

i ' : 

equations: ^ ' . . , . . 


Continuity : ~ ~ ~ ^o^^2 


;Cl3a) 


d>* + - (|> - 2(^=-^) 

^o ri o m 


^2 


^ - [qv ^ + 2 

0^0 1 o m 


*0 = ° 


Longitudinal momentum: (q - + n 


P ' 

11 


+ '’<^>('’2 -h 

t ^o 

i • 


npf + 2 (^^)f 

o m o 


1<2 - 


(13b) 


+ V 


nv' + 2(— ) v„ 
o m o 


= 0 



Lateral momentum: 




- 3(^) 
m 


<^2 + 


(n - 

O c. 


4 . . ^ ) (p 


^2 . 
r ^ 72 ^2 


V = 0 

o 


Energy : (iji ~ 4>^) 


?// p ^ ,2 

.0 


F’F^ 

F^ 




+ 2(^)(Y 72 ^ ^ ( !o ^ ^ 

I m «o ^o *0 ■ 


t T P* ■ 

2(izl!*) _ n(Y ^ - -^) 
^ m > '"'Tf to Fo 


Vo = 0 


(13c) 


KlSd) 


i u 

Except for additional terms corresponding to the order-o 

} 

terms retained in the normalized equations (3)^ these equa- 

! ='■ ' ■ ■ . ■ ■ ■ 

ti-ons- are very similar to the order-e equations; many of 

the coefficients are the same, and the only body shape 

I 

parameter that appears is m. The similarity functions 

j 

I 

V 2 (n), and 4 > 2 (n)> which form the solutions to 

these equations, will therefore be families of universal 

functions in the same sense as the other solutions are. 

Furthermore, Just as the order— e equations are independent 

2 

of the body slenderness perturbation (in 6 ), these equa- 

t 

tions are independent of the shock wave strength perturbation 
1 

j(in e). Thus, application of these equations to determine 

j 

the body slenderness perturbation of the zeroth-order small 

} ■ 

disturbance similar solutions neither requires nor excludes 
( 

application of the equations for. the shock strength 



21 


perturbation at the same time. Figure 2(c) shows that v/ith 

both perturbations applied the expected error for 

a given Mach number and body slenderness is much less than 
; 

without them (Figure 2(a)) or with just the shock strength 

perturbation (Figure 2(b)). 
i • 

\ Since the order-e solutions have been found previously 

( 

I O 

and are not needed to get the order~6 solutions, they V7ill 

I 

hot be considered further. All subsequent development will 

O O O ll 

assume e < 6 , so that e < e6 <6 << 1; all terms of 

i 4 

order 6 or smaller will be neglected. 

i 

i . ■ ' . 

j C . Boundary Conditions 

! 

j This section will deal with the boundary conditions at 
the shock and at the body surface and with the implications 
of the body boundary conditions on the solutions near the 

I 

surface. 

( 

; Shock wave . The boundary conditions at the shock wave 
are determined by the oblique shock relations (equation (4)). 
Using the expression for R(-^) (equation (7)) these rela- 
tions become: 


! 
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2 ^2 ^-2(l-m) . ,2 2 


P - "~ 7 T C 

Y+1 ^ 


|2(M:.2) _ i] 5 


-ij (l-rn) 


a, + 0('S‘*) 


= m L e _2_ (1+1) ^2(l^m) + 

Y-1 Y+1 Y-1^ -^2 ^ ■ 


S Y- 


u = — tt ^ C 
i s Y+1 


2 „2 _-2(l-m) .2 2 


r2(3lL_2) a - ijm** 5 
L m 2 J 


4 ,-4(l-m) 




2 ^ ^-(l-m) . ^2 2 r,3m-2, ^ ,1^3 r 

; s Y+1 Y+1 L m 2 J 


-3(l-m) 


; 2 r/2-m\ „ IT ^1-tn J. ncr*t\ 

+ ® ^ - sr ^ ^ 

' ^ m J 


bomparlng these equations term-by-term with equations (8) 
determines the boundary conditions for the similarity func- 
tions at the shock wave (n = n ): 

s 


F (n ) 
o^ s' 


^ (h ) 
^o s 














Y+1 2 1 

Y-1 Y-1 „2 
' m 

- -|t[ 2(— )a, - 
Y+1 m l 2-’ 

m 

2 r ^ 2 — m 1 1 

- -2^ 

' m 
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'*'2^'^s^ = 0 


''2("s> = - 7frf2(3^)a2 - 1] 


■f2(ng)' 


2 rf 3m-2 N _ ni 

Y+1^^ m ' 2 


\ 

>• (15c) 


J 


Note that the shock wave ciisplaceinent constants a^ and 3.2 

are Initially unknov/n. They depend on the parameter m and 

! 

are to be found in satisfying the boundary condition at the 

1 

body surface as part of the solution of the flow equations. 


! Body surface. The boundary conditions at the body 

\ ' 

} ■ ' 

surface are determined by the mass flow through the surface. 

If V is the velocity and p is the density of the flow 
w w 

put through the surface (Figure 3) y the mass flow balance 
hormal to the surface is given by: 


Pb% ®b - Pb“b ®b = Pw'^w 


( 16 ) 


Or, in the normalized variables: 


v^ = a v^ sec 0^ + (1 + ^ (17) 


i 



Figure 3. Vector' diagram of flow at the body surface 


where « = . Now 

I oo . 


and 


tan 0. = 

b dx b 


=[ 


sec 0^ = 1 + 6 riv. x 


2^2 ^2 ,^-2(l-m)' 

_ 


1/2 




) ( 18 ) 


, ^ 1 .2 2 2 -2(l~m) . 

= 1+ 2 <5n^ni X ' + 0(o). 


J 


Putting these into equation (1?) along with the expressions 

i 

in the similarity functions for u and v from equations 
(8) gives the relation 
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- % 


mC - OLV 


w 


+ 6 


2\ r 


i ^ 

I I 

+ e + 0(6'^) = 


3.-3(l-™) _. „^ ^^2 -2Cl-m)i 

W d 


(19) 


if the flow through the surface has the particular form 
-d-m) 


av„ = 

j w ^w ^ 


5 where ■ is a constant , the boundary 
w 


conditions at the body surface are (from equation (19)) 


= °* 




( 20 ) 


<(>i(nb) = 0 


While the development of these boundary condition shows that 

mass flow through the surface can be accommodated without 
[ 

difficulty, the rest of this dissertation will be restricted 

to the no mass flow conditions, <j> =0. The resulting 

w 

boundary conditions are 
} ^ 


i 

( 



\ 

■l-gCnb) ' = 0. (21) 

i I , 

1 I 

I : ■' 

i 

I Initial magnitude checks . These boundary conditions 

can be used with the flow equations to provide some initial 

i 

checks of the order of magnitude of the similarity functions 

) 

) 

As stated in Section II-B, these functions have been 

! 

assumed to be of order unity. From the boundary conditions 

at the shock (equations (15)), this assumption appears 

! . • . • * . 

I 

justified there, for y not too close to one and if ^2 is 

hot too large, except that Having a function 

1 

become much less than one does not invalidate the procedure 

used in getting the equations so long as the function does 

\ 

not appear as the denominator of terms that one dropped as 

ii 

being negligibly small, i.e. of order 6 . Neither nor 

' 2 

any other order - 6 function is in the denominator of any 
term that is dropped. 

I , Solutions to the zeroth-ordef equations given by 

i 

Kubota (9) indicate that and remain of order one 

5 

from the shock to the body but that goes to zero at 

I - ° 

,the body surface when m < 1. This result could affect the 
validity of the solutions in the region where 


is small 



since 


does appear in the denominator of a number of 
terms of equations (10) - (13)- One of these is the zeroth- ^ 
order lateral momentum equation (11). Applying the boundary 
condition ^ this equation becomes . 




'b’ 


from v/hich F’ (ni_) = 0. The zeroth-order longitudinal momen- 
[ o D 

turn equation (11) multiplied through by is 


j (n - «|> ) ^ V ' + nP ’ + 2(^^^)(v + P^) “ 0, 

I ' o o m 0^0 o 

i 

i 

! ■ 

Using 4>o^^b) ” ^ this becomes 

i ° 

Since il^ (Uv.) = 0 and F (n^.) f 0, this requires v “ 
as U Thus there’ is a (non-physical) - singularity^ in the 

similarity solutions at the surface of the body h = . 

i One possible way to avoid the singularity at the body 

i 

surface is to reformulate the problem. The fact that 
which remains finite as n is the zeroth-order 

'similarity form of the longitudinal momentum suggests that 


* Kubota (9) and Mlrels (10, l6) do not encounter this problem 
because they omit the variable entirely. It only occurs 

(..in the longitudinal momentum equation, which they do not 
use, prefering a Bernoulli equation for the velocity. 



using momentum components, instead t)f ve3.ocity components 
fundamental variables might remove the singularity. This 
approach is pursued in the next section. 


as 


D. Momentum Variable Formulation 


I As was noted in the previous section, the zeroth-order 
i 

longitudinal velocity similarity function v is singular 
i ^ 

I 

at the body surface. However, the product ip v remains 

I 

finite, which suggests that a reformulation of the problem 
in terms of new variables might remove the singularity and 


allow the numerical Integration to proceed all the way to 


the body surface. The variables chosen for the reformulation 

t 

are the longitudinal momentum pu (which is expected to 

■1 

1 

behave like and so remain finite at the body surface), 

the lateral momentum pv, the pressure p, and the density 
p. In terms of these variables, the Inviscid flow equations 
(1) become: 


I 

bontinuity: dm . 3Fv . ^ m « ' (23a) 


Longitudinal momentum: 




pu 


Ip 

dx 


) + 


PV(5 


— apx . -2 3p ^ 


i(23b) 
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1 

Lateral momentum: 


pu(p ^ 


(23c) 


Energy : (pu + pv |p) (^) = 0 


‘(23d) 


I Normalization . These equations are normalized using 
I 

the expressions: 




X 

^ I 


P = > (P^) == 


^ ' a ^ P 


P P< 


pu 




p u 
^ 00 00 


6 (p u ) 

' ^ nrs rv\ ^ 


(pv) = 


) (2^) 


pv 


<5P u 


J 


!(Note that these expressions are the same as equations (2) 

t 

with the exception that (pu) here is the same as the previous 
p(l+6^u).) The normalized forms of the equations (23) are 

_f 

then: 


Conti: 


3(pu) . 9(pv) . ^ (pv) _ 
+ “3F“ ^ “7 


0 


(25a) 


Lonp;ltudinal momentum : (pu) jp - (pu) -|^j 


5(pu) 

dr 


+ (pv)[p - (pu 


. .2 2 3 P 

) aFj^ ^ P 3 i = 


(25b) 


0 
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Lateral momentum 


(pu) 


9 ( p V ) 


- (pv) 


iP 

9x 


(25c) 


* (pv) p - (PV) |£ 


p »= I? . 0 




Ene_r^;^: |(pu) ^ + (pv) ^ 


(^) = 0 


(25d) 


-» P 


Y 


I Similarity variables. Just as in the first section of 

! 

this chapter, the similarity forms of the momentum variables 

I 

are chosen using the relations for flow through an oblique 


shock wave as guides. Combining the density and velocity 
relations used previously (equations (^)) gives the follow- 

j 

ing relations for the normalized variables just behind an 
bblique shock: - - ' v • - . .. . 


;( pu) 


(pu) 


p u 

*^OC» CO 


(Y+1) sin^ 

(y-1) sln^ ~ ^ 


1 - 


2(M^ sin^ 0-1 
«> s 


(I^ + i^) + (^^ + 


m 


- 1) 


<=¥■ 


m: 


? p ? 

00 


)(1 + 6 ^ R '^) 


(Y+1) 


(26a) 


= l±i - 

Y-1 


Y-1 


R’^ + e 


g,(-Y ± ll.. (Rt ^ 0(6^) 
(y-1)^ 
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( pv) 


iR 


(pu). 


2 (y+ 1) n ^ sin^ 0 (M ^ sin^ 0^ -- 1) 

' CO S S 


(Y~1)M 

L 


? P 

sin 0 +2 

‘ s 


(Y+1)M^ 

J 


tan 0 


,2 


k ^6 


i — (1 + 6 ^ R '^)1 R ’ 


2.2 


(26b) 


[Xl ri 2 + 1 — (1 + ^2j^,2^1 ^ 6^R»^) 

L ^ $^]yi ^ J 


Y-1 


Rt „ ^2 j^r3 ^ ^ (R«)“^ + 0(6^) 


Y-1 


(y-1) 


The similarity forms of the momentum variables are taken to 
be : 


pu 


= (r) + 6^ 02(r)R'^ + e u^(r)(R’) 




-2 


and 


) ( 27 ) 


2 


pv = y^(r) R^ + 6 "^ ^ Vi 3 ^(r)(R*) 




Using the same perturbed, power-law shock wave shape as 
before (equation (7)), these become (for e<(S^) 


pu = UgCn) + 


+ e o,(n) ^ + 0(6^) 

m 




and 


pv 


= Uq(t)) + e y^(ri) i m)>(28) 


+ 0(6 ) 
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The similarity forms of p and p are unchang6^d from the 

previous formulation (equations (8)), Putting the similarity 

forms of these variables into the flow equations (25)aand 

using the chain rule for partial derivatives (equations (9)) 

as before!, produces: 

! ! 


Continuity : 


(p' - qo' + — p ) ^ 

O O X] o 


m 


+ 6^[p^ - qu^ - ’Jo + ^ ^ 

l-In^ .. . a .. 1 1 ^(l-2m) . 


(2 9a) 


I + e[p-! ~ qu* + 2(^) u, + ^ p J i C 


1 *‘^1 ‘ m " "'1 ' q ‘"1^ m ^ ^ 


= 0 


Longitudinal momentum : ” ^o^o^ 


-1 


+ <S {(nu^ - JJo)(^o’^2 “ 


i Vj2 _ (no^ - p^)o<] ^2 + 


+ [(nu^ - + n(v; - - 2 (^^).^^u^]U2 

- I'o + 2(^) P^]} 


(29b) 


+ e{(no^ - - Vf - t2(^) u2 


+ (nu^ - P„)u-] 


+ [(nv)^ - 


+ 2(-~) u^} i 2m) ^ ^ Q_ 


t 
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Lateral momentum; 


[(nu - ]i )(y - ip uM + - (i:— 51) || u ] 

o o ^ 0^0 o o o m o o 


+ «^{(no^ - ■" ’‘’ 0^2 


[2i|j p» - (nu “ y )y^ + 

^00 o 0 0 m o o ^2 


1-m 


;+ [(nu^ - - 3(-^) Vo’>^2 


1-m 


+ [n(y ^ y’) - y ]u<^}m ^ 

1 ^^^00 00 m ^00 2 ^ 


4 ^3m-4 


!+ e{(nu^ - yoXu^iCi - ’I'oM].) + V-oPi 


+ [2^0 - 


+ [(^y^ “ y ) \|/' + (\|/ y’ - y iJjM + u ]y^ 

'• ' ' o 0 0 ^0 0 0^0 m o 1 


^ [nCUo-K; - v;) - + 0(6**) = 0 


I 


(29c) 

A 



Energy : [(nu 


o 


y ) (yF ij; ’ 
O 0^0 


i/j P' ) 
o o 


2(i- ~)P u ] 

in o^'o o ^ 


+ 6"^{(nu^ - [y(nu^ 


y )^’ 


1-m 


;- *t(^) F2 ■ M 


o’^o o' *^2 


r [(nu^ - w„) f; - 2(^)(y-i)f^uJ 


:+ [n(YF^^^ - 4,^P^) - 2(1^) u^} m5c'*'"-5 


X29(i) 


e{(no^ - u^)(yF^i|)| - i|/qF{) + Y(nu^ - Uo)<I'qFj^ 


- Yi^-)F- 4 2(1^) (y+1)F^0^]^^ 


1-m 


;+ [n(YP Ip’ - iji P') - F Ip ]u-, 

' 0^0 ^00 ' m ^ 0^0 1 


- (yFqI(<^ - <(-^F^) mr^ + 0(5'*) = 0 


Zeroth-order equations . Using the same reasoning as was 

employed in Section II-B^ the coefficients of the zeroth- 

2 

order term, the order-6 term, and the order-e term of each 
of these four equations must be equal to zero. The 
zeroth-order system of equations which results from 
recognizing this fact is: 
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Contin uity 


)J’ - riu* + ~ M =-• 0 
o o n o 


Longitudinal momentum: (uu - = 0 

a V Q o o o o o 


Lateral momentum 


( 30 ) 


,(nu_ - y^)(u ^ y’) + F’ 

> o o o o o o o o 


ij; u y = 0 
m ^o o o 


Energy: (yu - y )(yP - t]; P’ ) ~ 2(^^— ^)P u = 0 

o ^o ’ 0^0 o m 0^0 o 




The longitudinal momentum equation can be integrated 
Immediately. For (yu^ - y^) ^ 0 (i.e. away from the body 
surface; see next section), it becomes "* 

which has the solution, Comparing equations (26) 

and (28), the boundary condition on at the shock wave 

is ^^(n^) = Thus, with (n ) given by equation (l6) , 

the constant c is 



c 

u 

o 


*7^ 


(l±i) 

_ Y-l'^ 

(l±i) 



(31) 


(This result could have been anticipated by making a 
comparison of the similarity variables for the two 
formulations as given by equations (8) and (28).) 

Using the results just obtained, the three remaining 
zeroth-order equations may be written: 
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Continuity : 
Lateral momentum: 


w; - % = 0 


•(n^ - y )(]i U; ' - iJj y ’ ) + [F’ - ( ^ - ^ ) y ] = 0 


Ensi£l- - MoXyPo'I'; - ’•’o^o = ° 


> (32) 


J 


\ Order-6'^ equations . Also using the result ~ 

■ 2 

the order-6 system of equations may be written: 


Continuity : y^ - - 2(^^) ^^2 n ^2 ” ^ 


(33a) 


'Longitudinal momentum: 


Cnl/o - - U„)'p^]'p2 


(33b) 


+ [(n^„ - w„)4.; - 2(^)4 -^]o 2 - ^2[nP; + 2(i^)F^] = 0 


1-m 


I 


' Lateral momentum : (nip^ - ^o^^^o^2 ** ^^0^2 


+ [t„(2F' + - (n^„ - Ko^>'P^2 


(33o) 


+ [(n4i^ - u„)<l>: + (t„u: - u„K/0 - 3(ijp)'('2]u2 


o o ^oo oo 


1-m 


+ [ri(y 1 /)^ - i/; y') - i(> y.]i>o “ 

o o o o m 0 2 


0 



ii3 


Enersx.: - Uo)(yFo'I'^ - V-oPp - ~ Vo^'^2 

+ [yCnV-o - - 'J(^)>l'o]P2 

; i'(33d) 

I 

I - [( nn -„ - v ^) p ; - 2 (^)( y - i ) p ^<^^]./-2 

i 1 

I + - ’*'0^0 ) - 2(^)Fo'^o’^ 2 = ° i . 

J 

|The order-G equations governing the shock wave strength 
•perturbation are very similar and could be written in the 
pame manner, but since they are not needed here they will 

I 

'not be derived. 

; Since the physical momentum component variables are 
'simply the products of the density with the velocity com- 
ponents in the usual formulation, the similarity functions 
of the tvjo formulations are simply related. Noting that 
the normalized momentum variable (pu) is equivalent to 
ip(l + 6 u) in the previous formulation, the relations between 
'the similarity functions are: 

I 

f 

I 

1% " K’ *"o = Vo’ '^2 " *2 - Vo’ ’'a " Vs * '*'o'*’2 

I O'!) 

’ ' Boundary co_nditions . To determine the boundary condi- 

.tions for the momentum variables, equations (28) must be 



■ . , 1 

1 

compared to the expanded form of equations (26), Putting 
the shock wave shape (equation (7)) irito equation (26) and 
expanding gives 


■(pu) = lU - 6^ ^ ^ 

' s -r-1 -Y-l (Y_1)2 nj2 

I 

! (35) 

I 

l(pv)3 = [a2(3^) - 

: + = ^ U^)a, + ^] i + 0(6h 


[The term-by-term comparison of these equations with equations 

I 

( 

!(28) gives 




U2<np 


2 

Y-*!’ 



) = 


2 

Y’l 




1 ] 

(36) 


The boundary conditions at the shock on ^q’ ^2 

■are the same as given by equation (15). The boundary condi- 
tions at the body are determined by the mass flow through 
;the surface, as in the previous formulation. Section II-C, 
•In this case, when equation (16) is normalized in terms of 
I the momentum variables it becomes; 


(pv),^ 


(pv) sec 
w 



tan 0^ 


(37) 



Putting in the similarity functions^ (equations (P8)), noting 
u = dj ^ and the relations for tan 6, and sec 0, (equa- 
tions (l8))a the mass flow balance becomes: 


I 

' r / X / x\ 1 ^(1-m) / X rn . r.2 1 2 2j.-2(l-m) 


;+ 0(6 ) = 0 , 

If the mass flow through the surface has the form 
— ( l—in ) 

(pv)^ - mC ' \ where is a constant (y^ = 0 

2 

for no flow), the zeroth-order and order-6 boundary 
conditions are 




43 ( 0 ^) 


(n,_) 

b^o b 






1 r.2 

2 


(38) 


For no mass flow through the body surface, which is the case 
of present interest, these boundary conditions become 


" '’b'^'o^'^b^ ' >^o(nK) - HhOoCOv,) = 0 


2' 'b' 'b 2' 'b 


( 39 ) 


The boundary conditions can be used with the flow 
equations as in the previous formulation to provide some 



initial checks of the order of magnitude of the similarity 
variables. At the shock vjave, the boundary condition.s 
(equations (35)) Indicate the variables are of order unity 
there. At the body, the zero-order boundary condition 
(equation (39)) produces the same results from the zero- 
jorder energy equation ( 32 ) as before; i.e. i^^(rij^) = 0. 

Thus the boundary condition at the body (equation (39)) 

'becomes Vt^(n^) = ^b^o^^b^ ^o^^b ^ ^ ^b 

agreement with equations (3^) and (21). 

1 

The development above shows that it should be possible 
■ 2 

to get order-6 similarity solutions using the momentum 
^variable formulation, and that this formulation avoids the 
'Singularity in the . zeroth-order longitudinal velocity at the 
body surface. Whether the formulation is successful in 
^voiding singularities in the order-6 variables must be 
determined from the numerical solutions to be obtained. 


E . Correlation of Solutions 
Hornung (26) suggested a method to nondimensionalize 
experimental data so that shock wave shapes and pressure 
distributions for a given body power-law exponent would each 

form a single correlation curve. The basis of his correla- 
tion is to nondimensionalize the shock coordinates by a 

length scale D defined such that the body shape is given 
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For = njjx'”, D = (SHjj) ^ -S’. 

' ! 2 . 
jWhen this correlation is applied to the order-6 similarity 

isolation shock-wave shape, given by 



■Since this equation does not contain the slenderness parameter 

‘6, it gives a single curve for any given value of the body 

2 . 

power-law exponent m; i.e., the order-6 similarity 
‘solution produces a single correlation curve independent 
of the value of the slenderness parameter 6. 

Hornung suggested that pressure data would correlate in 
-much the same way if it were plotted as P/P^jo^^oo x/D, 

.where D is the same non-dimensionalizing scale length 

i 

2 

as used for the shocks. In fact, the order-6 solutions 
•show that not only the pressure, but also the other flow 
'Variables correlate when x/D is the longitudinal distance. 
'The pressure, density and velocity components from 
'equations (8) can be written 





^8 


P_ = 2(5^P ^ 26^ 


*n ^ \ \. i" TTl ' Z' \ 

F^m {j) 


+ ‘5^F2rn (|) 


= 2 


? P — P X “■2(1’’^) Ij h -4 X 

F^m^6'^(6ri^) (p) + ^ 


2 2(l-m) 

2F (5^) (i) 


^ 2F2(^) (|) 


4 - -4(l-ni) 


P I 4- ^ 2 /Xv 

^ = p = + 6 ii.2m (|;) 


2 ~ -*-2(l-*m) 
b 


^ = 1 + 6^u =1+6^ 

<x> 


P - -2(l-m) p y - -')(l-m) 
Vom2(|) + a^v^tn^l) 


‘2 - -2(l-m) 

= 1 ^ 


. . /tn ., /X V 


4 -4(l-m) 


and 


^ = 6v = 6 
u„ 


- -d-ra) P - -3(l-m) 

4'o™(|) + 6‘^4.2m^(|) 


- -d-ra) 


3 - -3d-m) 


= (P> 




[Since none of these contains 6 explicitly, each one forms 
a single curve for a given value of the body power-law 
exponent m. Independent of the body slenderness. Thus it 

i 

•should be possible to correlate experimental velocity 
component and density distributions, as v;ell as the pressure 
distributions, by plotting them against the normalized 
coordinate x/D. 



The regularity of the correlation form for expressing 
the physical variables in terras of the similarity functions 
suggests a possible refinement of the similarity formulation 
By using the local zeroth-order shock wave slope, expressed 
as 


I 

i 


I 




(l~m) 


:as the small parameter instead of the average shock wave 

f ' 

1 

■ ‘slope 6 = R^(I)/£ (Figure 1, p. 13) > it should be possible 
-to improve the formulation. In particular, this change 
Jwould facilitate the estimation of error in the nose region, 
'Where the local shock slope increases rapidly. While it has 
'not been possible to include it in the present study, such 
a reformulation would provide a good starting point for 

i 

further work on similarity solutions in hypersonic flow. 



CHAPTER III 


SOLUTION OP EQUATIONS 


A. General Scheme of Solution 
Either the velocity-variable or the momentum- variable 

I 

similarity form of the flow equation ((11) and (13) or (32) 
and (33)) with the boundary conditions at the shock wave 
(equations (15) or (35)) and at the body surface (equations 
(21) or (39)) is sufficient to completely determine a 
solution for the flow field. However, since there is no 
general analytic form for the solution in either formulation, 
it must be found numerically for each case (i.e. for each 
value of the power- law exponent m with o ,set equal to 
either zero or one) . The general scheme for obtaining the 
numerical solution is to begin at the shock U = Rg? where the 
boundary conditions are known, and to integrate the similarity 
functions numerically toward the body, which is known to be 
reached when the zeroth-order boundary condition is 
satisfied; i.e., U = when 

V (n) 

The derivatives of the similarity functions, used for this 
integration, are found by solving the flow equations for 
them algebraically. Thus, from the zeroth-order equations 
(11), the derivatives are 
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The derivatives of the order- 6 '^ functions^ from equations 
(13)5 are 

i 

^^(n) = j[(n-$„) ^ (3 - Y ^ i 

r o , — O -* i 


1/) F P^ P 

+ (ri — d) ) (-^ — ^ — ) ;=; — + 2( ) z ::: — + yqv ' 
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^000 o 
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+ [(n-4>o) (n<f'; + 


4 > ) + n 5 ^ + 2 (y+i)( 

m ^o P ' 

o 


¥>]| 


4- 


O 


^2 = '^o + [ 3 (^) - <^;]f 

[ 1 > F ' 

H(})^ + (^^)4> V { + -r^ 
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^o 
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[<=¥*'■ 

F' 
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i(4ld) 


Similarly, the derivatives of the similarity functions for 
the momentum variable formulation are, from equations (32) 


and (33): 
1-m 
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+ 2(ij^)(n - ^)(<l>2 + ^2^ + ■ n^’i “ ijr^] ^ 
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(1 - ijr^ - ^^>'<'2 + ^'^“2 ■" ‘"2^J 

^ 'f'^niUg - ’('2) + Y(^) 
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Two major difficulties must be overcome in order to 
apply the scheme of integrating either equations (40) and 
(41) or equations (42) and (43) from the shock to the body. 
One difficulty is the singularity at the body surface 
apparent from the fact that the denominators of some of the 

terms of these equations approach zero as the Independent 

’ i ^ 

variable r\ approaches the surface value It was over- 

come by using an approximate analytic solution, developed by 
Mlrels (10) and described In Section B of the Appendix to 
calculate the value of and the zeroth-order similarity 

functions in the region very near the body surface. For 
reasons explained in Appendix Section C, the second-order 
variables are calculated at the body by extrapolation of 
the order-6 similarity functions. The extrapolation techni 
ques used are described in the next section. The second of 
the two major difficulties is associated with the fact that 
the problem is a two-point boundary value problem. This 
difficulty is manifest in the need to choose initially the 
correct value of the shock wave displacement parameter a.^ 
(equation (7)) in order to satisfy the order-5 boundary 
condition at the body surface when h = at the end of 
the Integration, The steps taken to deal with this diffi- 
culty are described in Section C below. 

I • • ■ 

\ 

H 

^ as n from equation (21), 



2 

B . Extrapolation of Orderr<S Functions 
Two simple extrapolation techniques were used to carry 

p 

the order-<S similarity functions the short distance from 
the last computed point to the body surface. For the 
velocity-variable formulation the extrapolation used for 

i • . ■ 

each of the functions ^2’ ^2’ ^2 ^ cubic 

function of n passing through three computed points of the 

f 

function and having zero curvature at the body surface. 

The points used to define the curve were the last computed 
point and two previously computed points. The nximber of 
steps between the points was the next integer larger than 
the distance between the last point and the body divided by 

the last step size. . - 

For the momentum-variable formulation, the functlions 

i 

^2 and U2 were extrapolated linearly to the body surface 

2 

using the last computed point and slopes. The order-6 
stream function 62 and its derivative were also calculated 
at the last point, using the momentum-variable form of 
equation (A 9 ), and 82 was also extrapolated linearly. The 
values obtained were then used in equations ( 3 ^) and (A8) 
to calculate 

C. Methods for Determining the Constant a2 

Since the constant a2 is initially undetermined, the 

2 

value of the order-6 shock wave similarity ordinate 

I 

The results are insensitive to the technique (Section IV-B). 
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ti„ = 1 + unknown and cannot be used to 

S c 

begin the integration toward the body surface. Also, this i 

shock ordinate varies with the longitudinal distance so 

that its use would require a separate solution for every C , 

value. The use of the zeroth-order shock ordinate n = 1 

as the starting point for the integration avoids these two 

problems in determining the initial value of n„ but requires 

that the boundary conditions be transferred to n = 1 from 

2 

n “ the order-6 shock position, where they are known, 
s 

This transfer is made by using the Taylor series expansions 
of the similarity forms of the flow variables about the 
point n = in the same way as Kubota (9) and Mir els (10) 

S . ■ 

did for the order-e perturbation. The Taylor series expan- 
sion of a general function g(n) about is 

s 


g(n) = + 


dn 


n + 

s 



(n - n 

e — + 

* 


Applying this expansion to equations (8) (with terms of order 
e neglected) and evaluating at n = 1 give 


p(C,l) = F^(ng)m^C “ a2F^(ilg)3m^r^^^~'^^ 

|(^^a) 


+ 0 ( 6 ^) 


. 1 

p(5.1) = *o(ng) + - a2>)-;(n3)lm2r20-n') + o(6^) (44b) 
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u(^,l) = v^(ng)m^C - a 2 V^ (n^ ) 


4 --4 ( l~m) 


(44c) 


+ 0 ( 6 ) 


v(^,l) = j^>^(n^)mC (1-m) ^ - 


O' 's 




+ 0 ( 6 ^) 


Similarly applying the Taylor series expansion to equations 
( 28 ) gives i 


(pu)(5,l) = %(ng) + 6 ^[U 2 (Hs) - a 2 u;(ng)]m^C 2 (l-m) + oCa**) 


(45) 


(pv)(£;,l) = U„(n„)m? O-ra) + 6^[P2(tig) - a2VQ(ng)]m35 


o s 


+ 0(6 ) 


Evaluating equations ( 8 ) at n = 1 directly and comparing 
the results to equations (44) and (45) yields boundary 
condition transfer relations of the form: 


iI'qU) = 'l'o(ng), '(' 2 ( 1 ) = 'i'2(is^ “ ^2'*'o^''s^’ 


etc . 
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Putting equations (40) in for the derivatives and applying 

the boundary conditions at n (equations (l6)) gives the 

s 

transferred boundary conditions at n = 1- 


Fo(l) = - 7+1 

i 

>'2<i> - ^Ir pi]*2 - X 

♦2<1> - - ^ I3(3?)(S) - Ola, 


v^d) = 


2 )prf2ffi=2, _ _1_ + _2Xi„ - 1 

7+T m ^ ^ m ^ Y-1 Y+l' 2 


■tiod) = ^ |[(^^^) - 3(— ) + 2 -^la, - 1 
^2 Y+1 m ■ m Y+1 2 


Y+1" 2 


> (')6) 


Similarly for the momentum variables (equations (42) and(36) 






[3(^)(?^) - a]E2 + ll 


( 47 ) 


‘‘2(1) = ^ 2c]a2 - 1 
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Iteration method * These transferred boundary conditions 
provide a definite starting position for the integration 
toward the body, but the constant ^2 must still be deter- 
mined. There are two methods for determining a 2 » The more 
obvious one is to guess the value of a. 2 , integrate toward 
the body (using the method given in the previous section to 
reach the surface), test the order-6 boundary condition 
at the surface, and repeat using improved guesses until the 
surface boundary condition is satisfied closely enough. 

The improved gusses for this iteration method were made 
using the method of chords, a finite difference approximation 
to, the well known Newton-Raphson method. (Note that the amount 
by which the boundary condition is not satisfied corresponds 
to the mass flow through the surface according to equations 
(20) or (38)..) 

Decomposition method . The other method for determining 

a 2 takes advantage of the linearity of the equations in 

the order-6 functions, which allows superposition of 

solutions. It was used by Sakurai (4), Kubota (.9)j and 

Mirels (10) in obtaining their results and is applied in a 

2 

similar manner here. Each of the order-6 similarity 
functions is decomposed into a linear combination in the 
parameter a 2 ; 
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F^Cn) = ^ 


(^ 8 ) 

Splitting ■ each of the equations obtained in this way into 
two separate equations ^ by setting the term containing a 2 
and the other term each equal to zero, produces a system of 
equations in the subscript~a functions and a system in the 
subscript-c functions. The system in the subscript-c 
functions is identical to the original system of equations 
(13) or (33)* The system in the subscript-a functions is 
the same except that the Inhomogeneous terms (i.e. the terms 
that do not contain an order-6 function or its derivative) 
do not appear. These two systems of equations have two 
different sets of boundary conditions. In order to obtain 
them, the boundary conditions at n = 1 (equations (46) or 
(47)) are decomposed by comparisons with equation (48), 
giving 


2a^ ^ Y+1 m ^ m Y-1^ Y+1^ 2c'^ ^ 


Y+1 


(49a) 


- o], ^2^(1) = 


(49b) 




Y+1 


r ^ 3 , 41 - 2 ^) _ 1 ■■) + V (1) = 

m ^ ^ m ^W-1^ Y+l-*’ 


Y+1 


(49c) 
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y+1 


(49d) 


or 


= - v!t - a], U2,(l) = - 

^ - 3(^)(^^) + 2°J> ^2,(1) = - 


(50) 


2 

Y“1 


These decomposed variables are then substituted Into the 

O ' 

order-6 flow equations, so that the continuity equation (13)j 
for example, becomes: 


po^2a ^o n ^o'^2a ^o ^2a 

+ [-fi + ^ 4-0 + 2(^)H2aja2 

+ jvac (’''i + n ’^0^^2c - - ‘^0^*20 

+ t^o + ^ *0 + 2(i^)]f2c,- 'I Vi 
- [nv^ + 2(^)'>o]’^oj = ° 


Beginning at n = 1 with these boundary conditions, the 

1 

decomposed system of equations (in either the velocity 
variables or the momentum variables) is integrated toward 



the body. Near the surface the method given in Section III-B 
is used to obtain values for the decomposed functions at the 
surface. The boundary condition at the surface, expressed 
in terms j^f these surface values of the decomposed functions. 
Is (from equation (21)) 

I 

or (from equation (39)) 

[V2a(’1b^ - ~ V2c^\^’ = ° 

Thus, in the velocity-variable formulation, the value of a .2 
is found from the relation 


-2 


( 51 ) 


In the momentum formulation, is 


^2 = - 




%'^2a 


(iw) - P 




( 52 ) 


The value of ^2 is now used to recombine the decomposed 
similarity functions using relations such as equations (^8). 

Once these functions have been computed for any value 
of the body power-law exponent m (with a = 0 or 1), they 
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can be used to calculate the complete flovi field about 

any such body as long as it is slender enough that <S << 1 ; 

1 ' 

and the Mach number is large enough that e = 2 ~ 2 

I 

D, Description of the Numerical Method 
The equations der.ived in the previous sections of this 
dissertation have been programmed for numerical solution on 
the CDC 6600 digital computer at the NASA Langley Research 
Center. Three separate programs were written, correspond- 
ing to three of the different methods of obtaining solutions 
which have been discussed. Two of these programs integrate 
the velocity-function equations (^ 0 ) and (^ 1 ); one uses the 
iterative and the other the decomposition method for obtain- 
ing the. value of a.2* The third program integrates the 
momentum function equations ( 42 ) and ( 43 ) and determines a2 
by the iterative method. All three of these computer pro- 
grams use a standard Integration subroutine employing the 
fourth-order Runga-Kutta formula supplemented by a Richard- 
son’s extrapolation. This subroutine halves or doubles the 
integration step size automatically in order to meet a 
specified local truncation error. 

For the present computations the initial step size (in 
—ft 

n.) was 2 (.00390625), and the maximum allowable step size 

was 2”*^ (.OO7S125). Generally the step size decreased to 
-1ft 

less than 2 near the body. At each step estimates of 
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and were computed by tine method given in 

Appendix Section B (equations (Al8)), When both estimates 
agreed to within l.OxlO”^ on successive steps, the estimates 
were accepted as the actual values of and ^^b^ 

the values of the other functions at the body were computed 
from the approximate analytic solution given in Appendix 
Section B or the extrapolations given in Section The 

iterative programs used a "method of chords” algorithm to 
compute improved estimates of the values of a 2 . (This is 
a finite difference approximation to the well-known -Newton- 
Raphson method)* The iteration was considered to have con- 
verged when the order -6 boundary condition at the surface 
(equation (21) or (39)) was satisfied to within 0*5 x 10"*^^. 


CHAPTER IV. DISCUSSION OF RESULTS 


A. Zero - Order Functions 

The methods given in the last chapter have been used to. 

2 

compute the zeroth-order and order-6 similarity functions 

for a number of cases, which will be presented and discussed 

In this chapter.’ Unless otherwise noted, these cases are all 

for Y = 1.^5 representing air as an ideal gas. The Figures 

presenting the functions were plotted by Calcomp plotting 

machines directly from the computed results. The. slight 

wavlness which may be noticeable at some points in the 

Figures is a result of this computer-aided plotting process; 

however, the curves at all points on the plots are accurate 

to within 0.1 percent of the full scale values. 

The zeroth-order similarity functions F^, and 

are shown for several values of the power -law exponent 

m in Figure ^ for two-dimensional flow (a = 0) and Figure 5 

for axially symmetric flow (o = 1). These functions agree 

with the same functions calculated by Kubota (9) a Mlrels 

(10, l6), and Townsend (23). 

The pressure function F^ and the lateral-velocity 

function b are seen to be smooth and well-behaved from 
o 

the zero-order shock location (n = 1) to the body surface. 

i 

Note that the body surface values of lie on the line 

= n in accordance with the zeroth-order boundary condition 
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(b) Density function, 
Figure 4,*. Continued. 
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(a) Pressure function, P^Cn). 


Figure 5 \ 


Zeroth-order Similarity Functions for Axisym- 
raetric Flow (a = X). 
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(b) Density function, 
Figure 5. Continued, 
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(c) Longitudinal velocity function.^ Vq(ii). 
Figure 5. Continued, 
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(equation (21)). As expected from the discussion in Section 
II-C 5 the density function and the longitudinal velocity 

function exhibit different types of singular behavior at 

the body surface for m 7 ^ 1 ; goes to zero, and goes ^ 

toward minus infinity. From the Figure it appears that the 
effects of these singularities are confined to a thin layer 
near the body surface, at least for values of m near 1 . 

The zeroth-order lateral momentum function is 

shown in Figure 6 for both the planar and axisymmetric cases. 
As indicated in Section II^D, the longitudinal momentum . func- 
tion is the same as the density function and the 

pressure and density functions, F^ and are the same as 

in the velocity-variable formulation. In the separate cal- 
culations made for the momentum-variable formulations, these 
functions came out essentially identical to the previous 
calculations. Figures 4 and 5* The function is seen 

in Figure 6 to behave like 4^^, in agreement with the com- 
parison of methods in equations (34). Thus the reformulation 
in terms of momentum variables was successful' in avoiding 
the unlimited growth of near the body surface. 

B . Shock Displacement Constant 

2 

The variation of the calculated order -6 shock wave 
displacement constant a^ with the power-law exponent m 
is shown in Figure 7(a) for both two-dimensional and 



75 








77 


axlsynjjnetrlc bodies. The calculations made using the itera- 
tive method for obtaining 3^2 essentially the same 

values as were obtained by the decomposition method (equations 
(^8) - ( 52 )). As can be seen in the Figure, the momentum- ^ 
variable formulation results for a^ differ only slightly 
from the results of the velocity-variable formulation. 

However, all of the results are characterized by a singular 
discontinuity in 3^2 which has no physical counterpart in 
the actual flow about power-law bodies. This singularity, 
which occurs near m = .817 for a = 0 and near m = .653 
for a = 1, represents a rapid decrease in the distance 
from the shock to the body as the power-law decrease, 
followed by a Jump to a large distance at the discontinuity. 
Since this behavior is physically unrealistic, it must be an 
artifact of the mathematical processes (analytical and 
numerical) used to obtain the solutions. 

Figure 7(b) shows. that the singularity is associated 

with a zero in the denominator of equation (51)- (This is 

the equation used in the decomposition method to determine 

2 

the value of a^ satisfying the order-6 boundary condition 
at the surface.) Since there is a nearby zero in the 
numerators of this equation, it is probable that the 
singularity in the quotient a 2 occurs because the zeros 
in the numerator and denominator, which should coincide, are 
displaced relative to one another. (Coincidence of the 



Shock displocemenl cons^ont.Oz 


Figure 



Body power- low exponent, m 


(a) Computed values for a2. 

7 . Variation of Shock Displacement Constant b .2 with Body Power-Law 
Exponent m» 




(b) Numerator and denominator of expression for a 2 (equation (5D). 


Figure 7 ^ Concluded . 
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zeros would make a .2 mathematically indeterminate at that 
point but would allow a continuous variation of a .2 with m, 
from which the value at the Indeterminate point could be 


inferred.) The most likely causes of such a zero displace- 

! 

ment are accumulated truncation errors from the numerical 

I 


integration and errors in the extrapolation from the last 
iintegration point to the body surface. Any error occurring 
|ln these processes is aggravated by the fact that the 
[numerator of equation (51) is the small difference of two 


. Inearly equal extrapolated numbers. For example, with a 


0 


land m = .83> nv, == .73327 and the extrapolated values 
! ^ 
iVo(rjb^ = -9.9910 and 4^2c^^b^ “ -7.3153 give the numerator 

'-.0108 with- a relative error about 1000 times that of the 


■individual functions; that is, an 0.1^ error in 

\ 

;and ^2c^^b^ would become a 100^ error in the numerator of 

i ■ - / 

equation (51). 

; To test this line of reasoning, additional calculations 
[were made reducing the step size, reducing the range of the 


[extrapolation, and finally, extrapolating the whole numerator 

iof equations (51) rather than just the separate parts. The 
! 

Iresults of these calculations indicate that the value of ap 

I 

[is sensitive to these changes for power-law exponents less 

j 

jthan that at the singularity; but, the position of the 

I 

i singularity and the values of 3.2 . for m greater than that 

I 

tat the singularity were virtually unaffected. Thus, removing 
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the singularity would require a more radical change in the 
mathematical process than simply, changing parameters in the 
numerical integration and extrapolation schemes. This con- 
!clusion is supported by the results shown in Figure 7(a) 
where the singularity Itself appears to be unaffected by 

! 

« 

;the change to the momentum variable formulation and change 
•in extrapolation mode^ although there are large changes in 
;the computed value of a2 for m less than that at the 
^singularity^ Note that, while it is associated with the 

Izero in the denominator of equation (51) from the decom- 

i . ■ 

sposition method for determining a^a the singularity in 
|ap occurs at exactly the same place when the iteration 

t ■ ■ ■ ■■ ■■ 

method is used. 

i 

j The results of the present calculations should be good 

;for body power-law exponents above those for which there is 
i ■ • ' 

a significant influence from the singularity. Based on 

Figure 7, the range for good results is about .85 £ m £ 1.0 

for the two-dimensional flow (a = 0) and ,75 £ m < 1*0 for 

the axlsymmetric flow (a =. 1). Note that this is a more 

i 

restricted range than that for similarity solutions 
p 

! < m < 1, Lees and Kubota (8)) or to avoid first-order 

i — 

ientropy- layer effects ( ^ ^ 1 Freeman (17)). 
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2 

C . Qrder-6 Func»txons 
o 

The order-6 similarity functions ^2 

from the momentum-variable formulatioHj computed for axisym- 
metric flow (a = 1) using iteration, are shown in Figure 8. 

! I 

iThese functions are seen to have some differences in behavior 

Ifrom the zeroth-order functions described in Section IV-A. 

i 

♦ 

jThe most obvious difference is that the curves describing 
jthese functions do not all emanate from a single point at 

in = 1 . The variation at n = 1 is^due to the variation of 

1 

jthe boundary conditions at the shock with the body power-law 
|m and to the transfer of the boundary conditions from the 

shock to n = 1* As was noted Section B of Chapter III, the 

-I ' ' 2 

^values of these order-6 functions at the body surface were 

I 

i 

ifound using a linear extrapolation from the last computed 
ipoint of the numerical integration. The range of the 

I , ■ 

'extrapolation was less than .00002 in n for all the cases 
shown except for m = 1.0, for which the extrapolation was 
;made over an n distance of less than ,001. 

I 

i 2 

The order-6 similarity functions 4 ^ 2 » ^2* ^2 

ifrom the velocity variable formulation and using the cubic 
lextrapolatlon to obtain the surface values are shown in 
ipigure 9 for axisymmetric flow (a = 1). The results shown 

i 

iwere obtained using the decomposition technique for calculat- 

'ing the shock displacement constant a 2 (equations (48) - 

i 

1(52)), but essentially identical' results were also obtained 




(a) Pressure function, F2 (n)- 
2 

.Figure 8. Order-6 Similarity Functions from the Momentum 
i ' Variable Formulation for Axisymmetric Flow (a = 
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(b) Density function, 
Figure 8, Continued, 
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(c) Longitudinal momentum function, v) 2 (n). 
Figure 8, Continued, 
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(d) Lateral momentum function, 
Figure 8. Concluded. 
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(a) Pressure function, P 2 (n). 
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Figure 9. Order-6 Similarity Functions from the Velocity 
; • Variable Formulation for Axisymmetric Plow (a = 



5 


.80 


.85 


.90 


95 


(b) Density function, 
Figure 9. Continued. 
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(c) Longitudinal velocity function^ V 2 (n) . 
Figure 9. Continued, 
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(d) Lateral velocity function, cj^^Cn). 
. Figure 9. Concluded, 
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using iteration. The two functions and which are 

the same in both the velocity-variable and momentum-variably 
formulations, are practically the same in Figure 9 as in 
Figure 8. However, for power-law exponents less than those ^ 
■shown there are differences which become large as m de- 

I 

I 

jcreases, corresponding to the behavior of the constant a^ 

|( Figure 7). 

i The same comparison between methods of computing the body 

I 2 

•surface values of the order-6 similarity functions that was 

f 

made above for the axisymmetric bodies can be made for the 

•two-dimensional bodies (cr = 0) . The momentum-variable 

* 

i 

isimilarity functions for these bodies computed using the 
tiinear extrapolation technique are shown in Figure 10. The 

i 

jvelocity-variable similarity functions computed using the 

;cubic extrapolation technique are shown in Figure 11. These 

two sets of functions exhibit the same type of agreement as 

the functions for axisymmetric bodies (Figures 8 and 9)» 

: • 2 

All of the order-6 velocity and momentum functions 
;show singular behavior at the body surface for m < 1. It 

f 

was in anticipation of such behavior for that the 

momentum-variable formulation was undertaken (Section XI-D) . 

i 

1 

jWhile the reformulation was successful in avoiding the un- 

i 

;limited growth of at the body, it did not avoid the same 

1 2 

,type of behavior by the order-6 functions. Because of this 

i 

Ifailure to eliminate the singularities and because it does 




(a) Pressure function, P 2 (n). 


? 

Figure 10. Order-6 similarity functions from the velocity 
variable formulation for two-dimensional flow 
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(b) Density function, 

Figure 10. Continued. 
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(c) Longitudinal velocity function, v^Cn). 
Figure 10,. Continued, 
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(d) Lateral velocity function, <j> 2 (n). 
Figure 10. Concluded,. 





(a) Pressure function, Fp(n). 

' •••' . 2 

Figure 11. Order 6 similarity functions from the momentum' 
•variable formulation for two-dimensional flow 
(c = 0 ). 
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(c) Longitudinal momentum function, u^Cn). 
Figure 11. Continued. 


99 



(d) Lateral momentum function, p 2 (n). 
Figure 11. Concluded. 


o. 
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2 

not improve on the results of the order“*5 velocity-variable 
formulation in either quality or range of application, the 
momentum-variable formulation of the flow equations has very 
little advantage over the more conventional velocity-variable 
formulation. It has, however, provided a useful check of 
the numerical results. 

The singularity at the body surface is thought to be an 
entropy layer effect caused by the blunt nose of the body 
for m < 1. The effect is confined to a very narrow region 
since the very high curvature in the nose area (Infinite at 
X = 0) reduces the body slope rapidly. For example, when 
m = .80 and 6 = .5> the slope decreases from infinite to 
less than 1.0 before x/I = .006. 

As noted previously, Mirels (10, l6) recognised the 

singularity at the body surface and avoided calculating the 

longitudinal velocity function, which exhibits the singular 

behavior. In addition, Mirels (10) developed an approximate 

analytic solution which- is valid at the surface and used it 

to obtain surface values of the zeroth-order functions (see 

Appendix, Section A). The attempt to extend the approximate 

2 

analytic solution to order 6 is also discussed in the 
Appendix. 

Van Dyke (29, p. 185) notes the same singular behavior 
in a small disturbance solution at the surface for hypersonic 
flow over a blunted wedge and observes that it occurs because 



the zero-order solution '‘Is not a valid first approximation 
in the entropy layer.” He uses the method of matched asymp- 
totic expansions to obtain a uniformly valid analytic solu- 
tion for the blunted wedge. Adaptation of that method to 
the similarity solution problem for power-law bodies 
probably would extend the solutions to the body surface; but 
the application is not straightforward since the similar 
solutions are not analytic and therefore do not give the 
order of the singularities which occur. 

The effect of the singularity at the body surface, 

2 

then, is to make the order-6 solution inapplicable at the 
body surface except for the particular case m = 1, for 
which the body does not have a blunt nose. But, the 
singularity should not affect the solution away from the 
surface, where the similarity functions are of order unity, 
so long as the constant a .2 correctly determined. 

p 

Furthermore, the behavior of the order-6 similarity 
function for the pressure (P^) is quite regular all the 
way to the surface. Thus, the body surface pressure can 
be calculated using this function; however, the results 
must be suspect until checked against experiment or more 


exact results. 



102 


It is because of the singular ^behavior of the similarity 
functions that the numerical integration cannot proceed all 
the way to the body surface. Neither do the extrapolations 
follow the singular functions in giving values at the sur- 
face; so, I the calculated surface values of these functions 

i 

ho not represent the actual values of the singular functions, 

i 

which go to either plus or minus infinity at the surface. 

I 

.To the extent that they are useful for determining the value 
pf a2> the extrapolations can be considered as providing 


a limiting process for this constant. 

1 

! 


: Region of Validity of the Solutions 

I 

\ Three basic assumptions were required in order to obtain 

1 ' 

the hypersonic similarity solutions for power-law bodies: 

i 

;(1) the body is slender enough that terms of order 6^ are 

negligible compared to unity; (2) the shock wave about the 

■ 2 2 

body is strong enough that terms of order c = 1/6 M are 

j 

’negligible compared to unity; and (3) the Mach number is large 

'compared to unity. (The second of these can be relaxed to 
: 2 

e << 1 if the first-order solution in e is applied. ) 
However, even when these three assumptions are met overall 

|for a particular power-lav/ body, they generally are not all 

{ 

met In particular local regions, 
i ij ■ 

j The first assumption, 6 << 1, is obviously violated 

in the nose region of all blunt bodies, such as the power-law 
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bodies for ra < 1. Thus, the similarity solutions cannot be 

expected to apply at the nose of these bodies. But, the 
2 

order-6 solutions should be particularly useful in provid- 

I ' 

Ing an improved solution a moderate distance behind the nose. 
|As mentioned in the previous section, the violation of the 
{slender body assumption by the blunt nose is also the cause 

i . . • 

j ‘ 2 

jOf the singularities in the order-5 functions at the 
body surface. Therefore, because the first assumption is 

i 2 

^violated at the nose, the order-6 solution does not apply 

I 

at the body surface. 

I ’* 

I The strong shock assumption, e << 1> is violated wherever 

Ithe shock wave angle approaches the Mach angle, sin~*^(l/iyi ). 

^Unless the Mach number is extremely large, the shock wave 

♦ 

! 

will become weak far downstream from the nose of the body, 
land the similarity solutions will not apply in that region. 

I 

I The similarity solutions apply, then, in an Intermediate 

region from behind the nose to somewhere in the vicinity of 

' 2 

;the base of the body, and, in the case of the order-6 

i 

jsolution, only outside of the singularity at the body sur- 

! 

|face. The boundaries of this region depend on the Mach number 

j 

and on the power-law exponent and slenderness of the body. 

! 

'Decreasing the Mach number or increasing the slenderness of 

! 

I 

the body tends to weaken the shock; decreasing the power-law 
exponent increases the nose bluntness but weakens the shock 
•at the rear of the body. In any case, the boundaries of the 
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region In which the solutions apply are not sharply defined. 
They depend on the accuracy required in the results of each 
particular problem and must be ultimately determined by 
'comparison with experimental results, 

5 
? 

! 

I 

I The only exact solutions available for comparison with the 

‘Similarity solutions are those for flow over cones and 

i 

jwedges, corresponding to a power-law exponent of m = 1- 
i 2 

■Since the zeroth-order and order-6 similarity solutions do 

t 

!not contain any Mach number dependence, the most appropriate 

I 

jcomparison is at the hypersonic limit, ->• (To include 

i 

i 

iMach number effects the order-e solution would have to be 
iUsed also,) Results from the similarity solutions, in terms 

jof the physical flow variables, are compared in Figure 12 

! 

!wlth the exact solutions for flow over a range of wedge 
angles at Infinite Mach number. This flow, of course, is 
uniform behind the straight oblique shock wave, as indicated 

_by the solutions for m = 1 in Figure 4, pages 66 - 69? 

i 

'and Figures 10 and 11, pages 92-99. In Figure 13? results 

i 

'from the similarity solutions are compared with the exact 
^solutions for flow over circular cones.. Parts (a) of 

I 

(Figures 12 and 13 show the variation of the shock wave angle 

i 

(With the body surface angle, and include exact results from 
preferences 24 and 25 for ' several Mach numbers in addition to 


Comparison With Other Solutions 





(b) Pressure p/q^. 
Figure 12. Continued 





Magnitude of velocity .V — j — Angle between velocity ond 



Loterol velocity component, v/u« Longitudinal velocity component, u /u„ 



i (d) Velocity components. 

; Figure 12. Concluded. 



i , (a) Shock wave angle. 
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Figure 13. Comparison of Similarity Solution for ra = 1 
with Exact Solutions for Flow Over a Circular 
, Cone at jyi_ = 
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(c) Magnitude and direction of the velocity 
vector at the body surface. 

Figure 13, Continued. 
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Mc^ ” similarity results are found from 

_ , tan 6, 

tan 6 = T^t = 5R', where 6 = x == R’ is 

S 2n^,f Ht, 

found from equation (7) with m = 1,' The other parts of 

jthe figures show the variations vjith the body surface angle 

1 

pf the pressure coefficient p/Q^o? velocity 

- - - - - 

V/u^ = Vu + V /u^ , and the velocity components u/u^ 

and v/u in the uniform, flow behind the shock wave (Figure 

5 c» 

|l2) and at the body surface (Figure 13 ) • These similarity 

I ■ ' 

results are found from equations (8) with m = 1. The exact 
jresults are found from the oblique shock relations in the 

5 

f ■ • 

itwo~dimensional case and from the charts of reference 24 in 

;the axisymmetric case. 

I These Figures show that for m = 1 the zeroth-order 

! 

isimllarity solution agrees well with the exact solution for 

i 

body surface angles up to about 9^ = 12^, while the order*-6 

! ■ Q ' 

solution agrees well up to body angles of about 6^ = 20 . 

As can be seen at the bottom of the Figures, these cone or 

i - 

wedge angles correspond to slenderness parameter values of 
about 6 w .2 and 6 ~ .4. The similarity results for the 

f 

magnitude of the velocity show good agreement for even 

I 

I 

[larger body angles (Figures 12(c) and 13(c)). Since the 

i 

[error in the velocity components Is larger (Figures 12(d) 

i 

[and 13(d)), it must come largely from error in predicting 

j 

the direction of the velocity vector. This error in 

, direction is shown in the upper part of Figures 12(c) and 

• „ 2 4 

13 (c), where it is compared to the curves 6 9^ and -6 6^. 



These curves represent the order of error expected from 

2 

neglecting terms of order .6 and of order 6 , respectively. 
The error actually occurring is seen to be very close to 
that which was expected. It should be noted that this error 
jin the direction of the velocity vector cox-’responds to an 

j 

error in satisfying the boundary condition that there should 

I 

•be no flow through the body surface. That is, the velocity 

i 

Icomponent normal to the surface, ^ = y sin 

1 w 

should be zero (Figure 3, p. 30). Figures 12(c) and 13(c) 

2 

• jshow that this boundary condition is satisfied to order 6 
hy the order-6 similarity solution for m = 1. 

I In Parts (b) of Figures 12 and 13 > the pressure coeffi- 

L 

jcients are compared also with the Newtonian prediction: 

- _ 2 

•p/q^ = 2 sin 8^^. The Newtonian prediction is much more 

1 

laccurate for the conical flows (Figure 13(h)) than for the 

! 

jwedge flows (Figure 12(b)); however, even in the conical 
case the order-6 similarity solution is closer to the exact 
.-solution for moderately small values of the similarity 

i 

iparameter, i.e. 6£0.5* 

; In Figure 1^, the variation of the flow variables from 

the shock to the body is shown for three values of the 

! 

power-law exponent, m = -75, .85 and 1.0. The pressure, 

! 

■density and velocity components, calculated from the axisym- 

i 2 

■metric zeroth-order and order-6 similarity solutions, are 
‘shown at x/I = ,5 for a similarity parameter value of 


^ “l/bx . 
tan (— ) - e 

^b 
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6 = 0,4. The exact solution for a cone at Mach 20 ivith 
- 2 0^ (reference 25 ) is also shown in Figure l4(a) for 
comparison with the case m = 1.0. Although the conical 

bodies are not the same (0, = 20^ corresponds to 6 = . 3978 ), 

; i b . . 

2 

,the orderj-^ similarity results agree well with the Mach 20 
isolation shown for all the variables except the density. As 

j 

I 

[can be seen by the symbols representing the exact solution 
■for M ->' <» at the shock wave, the density is the only one 

;of the flow variables that is much affected by the difference 

] 

^between M = 20 and M -»■ ». The similarity solutions for 

' OO 00 

|the density agree exactly with the. infinite Mach number 
’ 2 

Isolation at the shock. The order-6 solution for the other 

I ■ . • ■ . : 

‘Variables differ from the exact , inf inlt e-Mach number solu- 

; ^4 

,tion at the shock by amounts which are of order o , as 

^expected from the approximation to the oblique shock relations 

used (equations (4)). On the other hand, the zeroth-order 

similarity solution is not accurate for a cone of this thlck- 

'ness; it is off by an amount of order 6 , which is l6^ for 

,6 = .4. Note that the zeroth-order results stop at q = 1, 

- 2 
which is the zeroth-order shock wave location. The order-6 

'results at the shock are given by equations (15) i the 

jintegrated solution begins at q = 1, with the results in 

1 

'the region 1 < n < n given by the Taylor series expansion 

s 

'as discussed in Section III-C, 
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Power-law exponent m - 0 . 85 * 
Figure l4 , Continued. 
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For body power-law exponents other than m = 
is no exact solution available for comparison with the order- 

p 

6 solution. However, there are some simple empirical method 
‘for estimating the pressure on general bodies . Two of these 

j 

methods will be used for comparisons. One is the Newtonian 
! 2 

ilaw Cp = 2 sin 6^. As discussed by Hayes and Probstein (1), 

j 

|it corresponds to the limits y 1 • 0 sind but, it 

\ 

t 

!ls widely used for more general hypersonic flows in this or 
modified form. The other empirical prediction is the 

jtangent-cone method, which takes as the pressure, at any point 

' i 

|on a body the pressure on the cone having the same surface 

i 

jangle as the body point. This method also is most accurate 
jf or 00 , since then the shock 'layer Is very • thin with 

jlittle pressure change across it. Hayes and Probstein (1) 

i 

■give a thorough discussion of these two methods and their 
ilimitations . Only one limitation will be mentioned here: 
these methods give only the body surface pressure and are 
;not complete flow field solutions, as are the similarity 

i 

isolations . 

f 

j 

; The calculated flow fields for bodies having m = .85 

( 

■and m = .75 are shown in Figures l^(b) and (c). In these 
I 2 

jcases the order-6 solution is again a major correction to 

i 

'the zeroth-order solution. However, the singularities 
; 2 

|in some of the order-6 similarity functions show up here 
’at the body surface. Because of the singularities, the 
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p 

order-6 values of the density and the velocity components 
are probably unrealistic close to the surface. Fortunately^ 
the pressure is well behaved all the way to the body surface, 
so that surface pressure coefficients can be calculated. 

! I ■ 

•; The surface pressures calculated from the similar 

j 

‘solutions are seen in Figures l^Ca), (b) and (c) to agree 

j 

ifairly well with the empirical predictions of the Newtonian 

i 

•and tangent-cone methods. However, as m goes from 1.0 to 
. 0.85 to 0.75 the order-6 similarity solution pressure 

i 

jdrops faster than the tangent-cone and Newtonian pressures, 

i 

\so that the agreement becomes progressively worse. Whether 
jthe similarity solution or the empirical methods give a 
tbetter representation of the actual pressure changes with body 
power law must be determined by comparison with experiment 
as in the following section. 

i ' ■ • ' 

f 

F, Comparison with Experimental Results 
: There is only a limited amount of useful experimental 

data on the hypersonic flow fields about power-law bodies. 

,'Thls data consists mainly of measured shock wave shapes and 

i 

surface pressure distributions for 3/^“ and 2/3-power bodies 
of a few different fineness ratios. In this section the 
, similarity solution predictions for shock wave shape and 
surface pressure for axisymmetric bodies with power-law 
exponent of m = 0.75 will be compared with the experimental 
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results of Kubota (9) and Peckham (19)* No comparison is 
made with results for smaller povjer-law exponents since 
valid order -6 similarity solutions were not obtained in 
those cases. Also, no comparisons are made with experimen- 
tal total drag measurements (e.g. those of reference 22 ) 

i 

because of the uncertainty in calculating the skin friction 

! 

contribution. 


. i Shock shape . Because the shock displacement constant 

!a 2 is so small, the order-6 shock wave shape is only 

.slightly different from the zeroth-order shape, an example 
i 

bf which is shown in Figure 1, page 13* For example, 
i 

bp = .0582 when m = .75 in axisymmetric flow (Figure 7 3 
page 78 , a = 1 ); thus, even for 6 = 0 .^, the order -6 

iterm in equation (7) for the shock shape R/£ amounts to 

f 

' only about .003 at the base of the body (x/I = 1,0). There- 
' 2 

fore the order -6 shock wave shape cannot be expected to be 


‘much closer to the experimentally measured shock shape than 

( 

|ls the zeroth-order shape. 

■ 2 

The zeroth-order and order -6 shock wave shape predic- 


|tions are compared to the shock shape data for m = .75 
■from references 9 and 19 in Figure 15* These data are pre- 

jsented in the correlation form used by Peckham (19) after a 

} 

'suggestion by Hornung (26). 



Dimensionless 



15* Shock wave shape correlation for axlsymmetrlc 3/^-power bodies. 
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The experimental data for the ,shock shapes are seen to 
correlate well with one another but to fall somewhat above 
the similarity solution predictions. This difference between 
experiment and theory is largely due to the difference in 

I I 

Mach numbers ( 6.85 and 7*7 for experiment, “ for 

1 ' 

■theory). It is comparable to the shift in shock location 
.With Mach number for cones (Figure 13(a), page 109). 

I 

I 

The experimental results correlate together because their 

I 

Mach numbers are relatively close. An additional cause for 
[the difference between the experiment and theory is the out- 

i 

Ward displacement of the flovj by the growth of the viscous 

boundary layer on the experimental bodies, 
i ■ ' 2 

P— —The effect of the order -6 term in the similarity 

j ■ • 

jsolution for the shock wave shape is seen in Figure 15 to 

'increase as x/D decreases. This is expected since small 

i _ _ _ _ 

Values of x/D correspond either to small values of x/^ 

or to large values of D/I, that is, either to points near 

the nose of the body where the slope is larger or to bodies 

which are less slender and thus have larger 6 values. Note 

J 

1 

■that the agreement of experiment with theory is better in 

I - - . 

'this region of smaller x/D values, as would be expected 

1 

i 

blnce the shock location is closer to that for infinite 

I 

Mach number for larger body slopes. The slope of Kubota’s 
■data for D/I = .0555 (circles in Figure 15) agrees very 

- 2 

'well with the slope of the order -6 similarity solution; 



this agreement In the slope on a log- log plot indicates good 
agreement of the power-lav/ exponent of the physical shock 
wave shape. 

• Additional shock wave shape data from hypersonic flows 

i - I 

over power-law bodies is presented by Freeman, Cash and 

t 

:Bedder (20) (m = .75s M ~ 8.8) and by Beaver (21) (m = »85> 
i ” 

^<x> " 7*0). Although they apparently correlate in the same 

I 

,way as the data in Figure 15 > these data are not presented 
*here because they fall in ranges of very large x/D values, 

i 

' 2 

!for which the order-6 term of the shock wave shape equation 

i 

'is negligible. For these large x/D values the strong shock 

j 

‘ 2 

■assumption, corresponding to e = 1/(M^6) << 1, may not be 

1 ' . . . , 

{satisfied. (See Section II-B). For example, e » .36 for 

I 

the 3/^-power bodies of reference 20. Thus the first-order 
■solution in e would have to be applied to obtain useful 
results. 

Pressure distribution . The pressure distributions 

obtained by Kubota (9) and Peckham (19) are shown in Figure 

16 for the same 3/^~power bodies as used for the shock 

wave shapes in Figure 15. .In addition to the zeroth-order 
2 

_and order-6 similarity solution predictions, the two 

I 

jempirical pressure distributions are presented. As seen in 
^Figure I 6 , the similarity solutions as well as the two 
^empirical methods give pressure distributions in good agree- 

j 

ment with the experimental data for the three bodies having 



Tongent -cone 
Zero-order similority 
Order- 8* similority 


Exper imeni 

O 8^0.555 f=l.03), 
□ 8= 0.269 f =2.13( 


,^r.7(ref 9) 


O 8=0.327, f-1.75,M„ = 6.85{ref 19) 


O 8=0.286, f = 2.00. M„ = 6.85 (ref 19) 


Comparison of similarity solutions with experimental 
pressure distributions; m = .75 (Theory is for same 
values of 6 as experiment; tangent-cone for same 
Mach number). 








For these cases 


fineness ratios f = l/2r^{i) of about two. 

2 

the order-6 ' similarity solution is very nearly the same as 
the zeroth-order solution except at the front of the body. 

The Newtonian prediction falls slightly higher than the 
:simllarity solutions back of the nose region, but curves 

representing these three methods are below the data points . 

! 

jsince at hypersonic speeds the viscous boundary layer 

I 

;tends to displace the flow outward, raising the pressure 

jabove that which would occur for Inviscid flow, the theore- 

\ 

J 

stical inviscid pressure levels are expected to fall slightly 

I 

^below those actually measured. For example, by applying a 

j 

iboundary layer displacement correction to the zeroth-order 
rsimilarity solution for the pressure on his fineness ratio 

i 

2.13 body, Kubota (9) obtained excellent agreement with his 
^experimental data (squares in Figure I 6 ). Since the tangent- 
cone pressure distribution falls slightly above the experi- 
mental data for the three higher-fineness-ratio bodies, the 
.other methods are somewhat preferable for these cases. 

' It is in the case of the body having a fineness ratio 
of nearly one, however, that a real difference between the 

methods appears. In particular, the difference between the 
’ 2 

, zeroth-order and order-6 pressure distributions becomes 
'substantial. The zeroth-order result lies above the 
lexperlmental data by an amount which Increases rapidly 
^toward the front of the body. On the other hand, the order-6 
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result lies below the data by an amount v/hich, for most of 
the body;* is only moderately larger than that for the finer 
bodies. This amount is on the order of the expected 

boundary layer displacement effect. It is only at the front 

; i 

s ! 2 

of the body that the order-6 result begins to diverge 
markedly from the experimental pressure distribution. Since 

the value of the slenderness parameter is 6 = 0.555 for 

\ ..... 2 • . 

this case* , it is not surprising that the order-6 similarl- 

i 

ty solution should begin to fail as the body surface angle 

1 

jincreases at the front of the body. This is about the same 

ivalue of 6 as the limit for good results in the wedge 
( 

land cone cases (Figures 12 and 13 > pages 105 - 112 ), The 
i . ■ 

‘[Newtonian method gives excellent agreement with the experi- 
mental data in this case, but this must be somewhat fortuitous 
'dn that no correction was made to account for the boundary 
layer displacement effect . The tangent-cone method again 
dies somewhat above the experimental data. 


*Kubota's (9) value 6 = .485 shown with his data corresponds 
i to 6n^ = .485 as used herein. 



CHAPTER V. CONCLUSIONS 


Beginning with the equations for conservation of mass, 
conservation of momentum, and conservation of energy for 
.the inviscid, two-dimensional or axisymmetric adiabatic flow^ 
.of an ideal gas, similarity solutions have been found which 

I 2 

give the flow field to order- 6 about povrer-law bodies in 
I'the hypersonic limit *> where 6 is a body slenderness 

l 

parameter* On the basis of this investigation the following 

1 

conclusions can be made: 

I 1. The order-6 solutions are independent of the 

^slenderness parameter 6. Thus the functions expressing 

i 

the solutions are universal in that they apply for all values 

L,. ■ ' ' 

’of 6 for which 6^ << 1. The relations between these 

! 

similarity functions and the physical flow variables are 
'relatively simple. 

2. For the purpose of obtaining the similarity solu- 
tions the flow equations can be formulated in terms of the 
•longitudinal and lateral momenta as basic variables instead 
of the corresponding velocity components, and the similarity 

1 

Ifunctions can be obtained. However, the expected advantage 
’of avoiding singular behavior of the equations at the body 

i 

■surface do^es not materialize. 

( . 

! 3 • In the present formulations the value of a. 2 > the 

^ 2 ■ 

shock displacement constant in the order-6 solution, goes 
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through plus and minus infinity at about m ~ .817 In the 

two-dimensional case and m = .653 in the axisyrnmetric 

case, where m is the body power-law exponent. Because the 

singularity does not correspond to actual flovj conditions 
i 

it must arise through the mathematical development. Since 
the singularity was not removed by any of the variations in 

I 

^solution procedure tried, the present results are limited to 
a range judged relatively free of effects from the singularity 
■(.85 £ m £ 1.0 for two-dimensional flow, .75 £ m £ for 

'axisyrnmetric flow) . 

; 4 . In comparisons with the exact solutions for inviscid 

2 

:flow over wedges and circular cones, the order-6 similarity 

} 

^results give excellent agreement for 5 less than about .4, 

corresponding to wedge or cone angles up to about 20^. Over 

2 

an even larger range, the order-6 surface pressure predic- 
tions were superior to the Newtonian pressure law. The 
2 

order-6 results were a significant improvement over the 

zeroth-order results for body angles greater than about 12*^. 

5* In comparisons with experimental shock wave shapes 

and surface pressure distributions for 3/4-power axisyrnmetric 

2 

bodies, the order-6 similarity solutions gave good results, 

considering that Mach number and boundary layer displacement 

effects are not included in the theory. For body fineness 

2 

ratios near two, the effects of the order-6 terms are 
significant only very near the body nose, whereas for a 
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fineness ratio near unity the order-6'^ terms had a large 
effect over almost the entire body. These good results for 
the surface pressure were obtained despite the singular 
behavior of some other variables at the surface. 

; 6. While the order-6 similarity solutions were 

I 

jdeveloped for the hypersonic limit the derivation 

jshows that they are compatible with the order-e solutions 
;of Kubota (9) and Mirels (l6), V7here e e 1/(M^6)^, The 

I 

brder-e solutions introduce Mach number effects* 

! 

I 7» While all present results were obtained for no 

) 

t 

jflow through the body surface as a boundary condition, it 
jappears from the derivation that small amounts of blowing 
W- suction through the wall could be easily accommodated. 

j 

j 

i 8 . It was noted that the correlation suggested by 

{ ' 

Hornung (26) for the shock wave shape and body pressure 

'distribution can be applied exactly to all of the flow 

2 

variables in the order-6 similarity solution form. This 
;finding suggests a method of correlating for future 
'experimental data. The form of correlation also suggests 
'a possible refinement of the derivation of the order-6*^ 
’similarity results, using the local body or shock wave slope 
:as the small parameter . 

I 

i 

1 

! 

I 

I 

I 

i 



BIBLIOGRAPHY 


1. Hayes, Wallace D.; and Probstein, Ronald P.: Hypersonic 

Flow Theory . (Second edition). Academic Press, 

New York, I 966 . 

2. Taylor, G, I. : The Formulation of a Blast Wave by a 

! Very Intense Explosion. Proc . Roy. Soc . London, 

! Ser. A, Vol. 201, 1950, pp. 159-186. 

3« Sakarai, Akira: On the Propagation and Structure of a 

; Blast Wave, I. J. Phys, Soc.. Japan, Vol. 8, no. 

Sept. - Oct. 1953 > PP • 662-669* 

A. Sakarai, Akira: On the Propagation and Structure of a 

‘ Blast Wave, II. J. Phys. Soc. Japan, Vol, 9 , no. 2, 

I March - April 195^ > PP* 256-266, 

[ 

‘ 5 . Hayes, W. D. : On Hypersonic Similitude. Quart . of 

: Appl. Math. Vol. 5, no. 1, April 19^7 > PP* 105-106. 

:6. Van Dyke, Milton, D.: A Study of Hypersonic Small- 

I Disturbance Theory. NACA Kept.' 119^ j 195^* 

7 . Lees, L. : Inviscid Hypersonic Flow over Blunt Nosed 

Slender Bodies. GALCIT Hypersonic Research Project, 
Memo No. Slj Feb. 1956. 

8. Lees, Lester; and Kubota, Toshi : Inviscid Hypersonic 

Flow over Blunt-Nosed Slender Bodies. J. Aero. > 

Sci., Vol. 2^, no. 3, March 1957, PP* 195-202, 

9 . Kubota, Toshi: Investigation of Flow around Simple 

; Bodies in Hypersonic Flow GALCIT Hypersonic Res. 

ProJ . Memo ^0, June 25, 1957* Also, Inviscid Hyper- 
: sonic Flow over Blunt-Nosed Slender Bodies. GALCIT 

' Publication No. 417, 1957* 

10. Mirels, Harold; Approximate Analytical Solutions for 

Hypersonic Flow over Slender Power Law Bodies, NASA 
! TR R- 15 , 1959 . 

11. Sedov, L. I.; On Certain Unsteady Motions of a 

. Compressible Fluid. Prlkl. Mat. Mekh. , Vol. 9, 1945, 

■ pp. 293-311' Transl. no. T-57 , The Rand Corp., 

Santa Monica, Calif., 1956, cited by Hayes and 
Probstein (1 ) , 


131 



132 


12. Sedov, L. I.: Propagation of Strong Blast VJaves. 

Prikl. Mat. Mekh., Vol, 10, 19^16, pp. 2U1-250, cited 
by Hayes and Probstein (1), . • . 

13 . Grodzovskii, G. L. : Certain Peculiarities of the 

Plow around Bodies at High Supersonic Velocities. Izv . 

, Akad Nauk SSSR, Otd. Tekhn. Nauk, no. 6, 1957s PP * , 

86~92 . (Transl . by Morris D. Friedman, Inc.) 

i 

1;^ . Chernyi, G. G.: Introduction to Hypersonic Flow . (R. F. 

^ Probstein, transl. and ed.) Academic Press, New York, 

; 1961 . 

i 

j 

15 . Stanyakovich , K. P. : Unsteady Motion of Continuous Media. 

f (J. G. Adashko, transl., and M. Holt, ed,), Pergamon 

M Press, New York, i 960 , 

1 6 . Mir els, Harold: Hypersonic Plow over Slender Power-Law 

I and Related Bodies. Advances in Applied Mechanics , 

j Vol. VII, Academic Press, New York, 1962, pp. 1-5^4 , 317 - 19 . 

i 

17 . Freeman, N. C. : Asymptotic Solutions in Hypersonic Plow; 

I An Approach to Second-Order Solutions of Hypersonic 

j Small Disturbance Theory. Research Frontiers in Fluid 

Dynamics (R. J . Seeger and G. Temple, eds.) Interscience 
I Monographs and Texts in Physics and Astronomy, Vol. XV, 

' John Wiley' and Sons, I 965 , PP • 284-307. 

18 . Sychev, V. V.: On the Theory of Hypersonic Plow over 

: Blunt-Nosed Slender Bodies. Advances in Aeronautical 

Science , vol. 3 (Th. von Karmdn, chmn. ed . comm . ) , 

' Macmillan Co., New York, 1962 , pp. 87-102. 

19 . Peckham, D. H. : Measurements of Pressure Distribution 

and Shock-Wave Shape on Power-Law Bodies at a Mach 
Number of 6 . 85 . Aeronautical Research Council C.P. 

871 , April 1965 . 

20. Freeman, N. C., Cash, R. F. , and Redder, D. : An 

J Experimental Investigation of Asymptotic Hypersonic 

Flows, J, Fluid Mech,, Vol. I 8 , no, 3, March 

! 1964, pp. 379-384. 

i 

21. Beavers, G. S. : Shock Wave Shapes on Hypersonic 

i Axisymmetric Power-Law Bodies. AIAA J., Vol. 7j 

* no. 10, Oct. 1969 , pp. 2038-2040. 

22. Spencer, Bernard, Jr,, and Fox, Charles H., Jr.: 

i Hypersonic Aerodynamic Performance of Minimum-Wave- 

Drag Bodies. NASA TR R-250, I 966 . 




23- Townsend, James C.: Hypersonic Aerodynamic Characteris- 

tics of a Family of Power -Law, V/ing-Body Conf ig*urat ions 
NASA TN D-7^27> 1973- 

24. Equations, Tables, and Charts for Compressible Flow. 

NACA Report 1135, 1953- 

25- Sims, Joseph L. : Tables of Supersonic Flow around 

Right, Circular Cones at Zero Angle of Attack. NASA 
SP-3004, 1964. 

26, Hornung, H. G.: Inviscid Flow over Plane Power Law 

Bodies. J. Fluid Mech., vol. 27, no. 2, Feb. 

.1967, pp. 315-336. 

27- Ashby, George C,, Jr.: Longitudinal Aerodynamic 

Performance of a Series of Pov/er-Law and Minimum 
Drag Bodies at Mach 6 and Several Reynolds Numbers. 

NASA TM X-2713, August 1974. 

28. Ashby, George C., Jr.; and Harris, Julius, E. : Boundary- 

Layer Transition and Displac.ement-Thickness Effects 
on Zero-Lift Drag of a Series of Power-Law Bodies at 
Mach 6,. NASA TN D-77 23 , ,, September 1974. 

29- Van Dyke, Milton: Perturbation Methods in Fluid 

Mechanics . Academic Press, Nevj York, 196^4". 

30. Van Dyke, Milton: Extension, Analysis, and Improvement 

of Perturbation Series, Tenth Naval Hydrodynamics 
Symp . , M.I.T., Cambridge, MA, June 27, 1974. 



APPENDIX. ASYMPTOTIC SOLUTION IN TERMS OF 


STREAM FUNCTIONS 


A . Stream Function Formulation 

i This! Appendix describes the development of an 

i I 

.asymptotically valid analytic solution to the flow equations 
^in terms of a similarity-function form of the stream func- 

j 

Ition. The solution is basically the zeroth-order solution 
; ? 

of Mirels (IO3I6) with order-6 terms added. It is 
developed from the velocity-variable formulation, which 
Mirels used, but the results can be easily related to the 

I 

•momentum-variable formulation also. 

i By definition, the stream function must satisfy the 
'continuity equation (3)5 which may be rewritten 

i 

i 

[p(l + 62 u)] + (pv) + a ^ =0. 


If 


and 


p(l + 6^u) 


1 ^ 

o dr 
r 


pv 


1 _ ^ 

3x 


(Al) 


/\ 

then i|j(x,r) is a stream function. It can be put into the 
s similarity form 
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0 ( n ) 

o 

1 + a 




(A 2 ) 


where the zeroth-order term was given by Mlrels ( 10 ) and 

2 

the constants B and J In the order~6 term are to be 


^determined by comparison with the similarity form of the 
iflov/ variables (equation (8)). Putting the stream function 

i 

1 

'defining relations (Al) into similarity form (using 

[ 

equation (A 2 ) and the chain rule equations ( 9 ))- 


p(l + 6^u) = + 6^B (l+a)m 


(l+a)n 


and 


(A 3 ) 


fp-v • =- ” 


[9 _ ^0.] Gl _ ^2 j-(l+ara) 

*■ o 1+a o o ■ m 2 2 a ^ 

n n 


The comparable relations from equations (8), omitting the 
'order-e terms, are 


p(l + 6^u) = to + '+ +0(6^) 

and • (AA) 


Comparing these two sets of equations, the exponent in the 

) 

1 P 

, order-6 term of (A2) must be J = (3 + <?)m - 2, and the 

j 

■similarity functions are related by the equations 
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o 


0^ 
o 

(l+o)n 


5 

o 


lll (j) 
^0^,0 


n 


n6’ 

_ (1^. 
a 1+0 




’Po + = 


o’ O 


B0^ 
m n 


V 2 + Va = -A tne^ ~ 


m n 


(A5) 


The boundary condition on the stream function similarity 
functions are 


e^(Pj,) = 0 , e„(i) = 1 


L,... 


e^Chb) = 0., e^d) = 1 


(A6) 


Solving equations (A5) for the constant B and evaluating 
'the similarity functions at n = 1, using equations (^6) 
and (A6), gives 

B = - ap 
y-1 2 


so that 


7cir £.(l+o)m a 2 ^^2 „ „2j,(3+a)m-2 

= 1+5- « - « ;pj e^m « 


(A7) 


:in addition, equations (A5) can be. solved for 
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:^2 = 


- 


( 1 + 0 )ti' 


K ° ^ w 


'^2 a 1+CT ®o ■'■ Y-1 ®'2®2J 

X] -• 


[ 

- Jf |»„»o * [< 


e 8 


(l+a) ^ - i3±«)lt2 e 
0' m 


(A8) 


I 


Or, conversely. 


= (n - 


O' ^O 






(A9) 


i As Mirels (10) showed for the zeroth-order case, the 
‘pressure can be related to the stream function by using the 
energy equation. Defining an entropy function to = p/p"'^, 
■the energy equation (3) in terms of the stream function 
becomes 


3ijj 8(0 _ ^ 

8r 3x 3x 9r 


(AlO) 


This has as a solution (to order 6"^) 


to = 


2_ (1= 


Y+1 Y+1 


e-® Lh, 

Y+1 o 


Y" °0 

(All) 


'] 
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where the zeroth-order term Is Hirels ' zeroth-order solu- 

% 

2 

tion (10) and the order-6 term was found by trial and 
error. The parameter 3 is defined by 3 = ~ 

2 

varies from zero to one as m varies from one to — , 

•the lower! limit for similarity solutions. The stream func- 

i - 

tion expression for the pressure can be found from the 
lentropy function, since the expression for the density is 

I 

•already given by equations (A8). Thus the density and pressure 
'are 


f Y 
;P 


- ,,.Y 


+ 


,2 ’<'2 
Y 


m) ^ 0(6^) ] 


r 9 ' IT 

= - — ^ r 1 - 

_(l+a)Ti^J _ 


^2 . . 2(l+a) ^ 


0' ^ 
ht) 


o J 


+ 0(6 ) 


and 


Y 

p = p'w 


_ 2 1 a 2 -2(l-m) 

- \±rr) (rr:^) ' o m C 


Y+1 Y+1 1+0 oy^.3 

n 0 

o 


(A12) 


I 


X <1 4- 6 


0 


^ - Y qt) - Trvjm^s 2(1 m)^ ^ ^^^4^ 


J „ 2 ^- 2 ( l - m )| 


|0n comparison with the similarity form of the normalized 
pressure (equation (8)), this equation gives, the relations 
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2 Y-l 1 
Y+1 '^Y+1^ 4+a^ ^oyjj 


n ' 0 


and 


o 


F 2 (n) = 33(6 ^ - Y ^) - Yv„(n)]Po(n) 


(A13) 


iwhere, again, the zeroth-order function is Mitels'. 

i 

i 

I Zeroth-Order Approximate Solution 

I In reference 10 , Mirels shows that, using equations (A8) 

and (AI 3 ) for > 4) 3 and P , the zeroth-order lateral 

^0 0 o 

'momentum equation (11) can be approximated for 8 (n) <<. 1 

i O 

hy 


6 " 

o 


3 


®o ^ ®o 


(Al4) 


He also gives the solution to this approximation: 


6„(n) = 


where 


- I’^r- 


(AI5) 


(A16) 


Mirels uses this solution to make an Improved approximation 

1 

i 

;for the lateral momentum equation and an improved approximate 
•solution 
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0 ( n ) 

. o 


Ko(n 


1+G 1+a 


__X_ 

Y-3 


- nv ) fl + 


l-a 

3 . 

2(2Y“3) 


(A17) 


Putting this second approximation back into equations (A8) 
'and (A13), he obtains the relations 


:n “ 


^ Y-3 


1+a 1+a 

n - n 


enj ° 6. 


'b ^"b “o i 

o ■ 2(2Y-S)F„(n^) ‘ 


& 


®o 


;^o ^ + r2Y-e ' ^,Cn^) -J 


(Al8) 


+ I ni ® e„ 

o O b 2 b o 


Using equations (A8) again, he obtains from these relations: 


n 


1+a 


= 1 - ( ^)[1 + 

Y-e ' n 


e(n-4'„) il< n. l-a 


2(2y-3) 


and 


(A19) 


F„(nv,) F (n) - I nv,(n-c(i )'(' (— ) 

o b o 2 b o ^o n, 

b 


These relations are correct to the order of 0 and 

o 

therefore apply near the body surface, where 0 (u, ) ™ 0 

i ^ ^ 

‘(equation (a 6)). Using (AIS), terms of order 0 are less 

y-3 

:■ 2 ^ 

than terms of order-6 when n-bv. < 

(l+o)n° 


. At the 


point where equations (AI 9 ) are used in the numerical 


solution to determine the values of q, and F (n, ), the 

b o b ^ 



integration has come so close to the body that 
Y-g 

(6/K ) ^ 

H-Hk — ‘ Thus the error in relations (A17) to 

(l+a)n^ 

*(A19) at Ithat point is much less than the order-6^ error 

I ! 

I i 

jof the equations being integrated, and the zeroth-order 
iapproximate solution for and (j)^ (relations (Al8)) 

:may be used the rest of the way to the body surface. 

• The zeroth“Order longitudinal momentum equation can be 

i 

iused to obtain an approximate solution for in the 

i . 

'region near the body. Noting the implication for 

Trom the body boundary conditions (equation (22)), a trial 

,} 

Isolution is made in the form 


Vo(n) 


v^(n) 


= D 


P^(n) 


(A20) 


F F' 
p - 9 . (- 9 . 

il; ^F 
^o o 


— ) 


.Putting this into the longitudinal momentum equation (11) 
gives 


D(n - 

O 


Ip ’ F 

iT^ + F 

^O O 


o + 2(^^)(D + 1) = 0. 

m 


Solving for D and substituting -from equations (A8) and 

! 

!(A13) : 



1^42 


3 + 


Y (''o 
1+a ^ 


0” 




6 ’ 

0 _o 

T 


n 


)n 


D = 


1 0” 

e(^) + (Y-i)(^ 

' O 


ft 0' 

R __0 
Y 0^ 


as ^ 


But D is a constant, so an approximation must be made; 

adopting Mirels first approximation (equation (Al4)) reduces 

p to “ that the relation (A20) for becomes 

‘(to order 9 , 6 << 1) 

, o ’ o 


v„(n) 


Y 

Y-1 


Fo(n) 

On! 


2 , 
Y -1 




(ria) 


Y 


(e^r 


(A21) 
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! ■ C. Order-6 Functions 

; In principle, expressions (A8), (A13) and (A21) could be 
‘ 2 

put ' into the order-6 lateral momentum equation ( 13 ) to 

■ 2 
obtain a first-order differential equation for the order-6 

similarity stream function However, the equation would 

contain a very large number of terms and would require a 

considerable amount of approximation to be made in order 

to reduce it to a form for which a solution could be 

obtained. The approximation process is complicated by the 

facts that 6^ occurs to non-integral powers and that many 

i 

terms contain 6, which Itself takes on very small values 

Xn the range of interest. Considering that the solution is 

only to be applied in a very small region near the body 

2 

surface and that Bg is already of order-6 , the additional 



accuracy vrhich might be obtained probably does not justify 

the additional effort required to produce an approximate 

solution in this way. Instead, the simple linear extrapola 

tion techniques described in Section III.B are used to 
• 2 

harry the order-6 similarity functions the short distance 
|to the body surface. 



